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THE  DIVERSITY  ECCM  PERFORMANCE  OF  FREQUENCY-HOPPING  CPFSK 
IN  PARTIAL-BAND  NOISE  JAMMING 

1.0  INTRODUCTION 
1.1  BACKGROUND 

While  frequency-hopped  MFSK  (M-ary  frequency-shift  keying)  waveforms  are 
widely  discussed  as  candidates  for  spread  spectrum  applications,  there  are  many 
current  applications  in  which  the  chosen  hop  modulation  is  narrowband  digital 
FM,  or  CPFSK  (continuous  phase  FSK)  in  its  many  forms,  most  notably  Army 
tactical  radios.  In  addition  to  offering  efficient  use  of  the  available 
spectrum,  it  has  been  estimated  that  limiter-discriminator  detection  of  a 
hopped  CPFSK  waveform  can  obtain  a  4  dB  improvement  in  performance  over  MFSK 
systems  in  noise  jamming  [21].  However,  the  few  published  analyses  giving 
results  for  FH/ CPFSK  are  either  based  on  approximations  or  neglect  thermal 
noise.  We  know  of  none  that  address  tone  jamming  of  FH/CPFSK. 

For  the  optimal  utilization  and  design  of  these  systems,  it  is  desirable 
to  discover  what  parameter  values  work  best  under  jamming.  Figure  1.1-1 
illustrates  a  slow-hopping  FH/CPFSK  system  without  diversity.  Shown  in  the 
figure  are  optional  transmitter  elements  associated  with  transmission  bandwidth 
control.  The  primary  transmission  parameter  is  the  normalized  maximum  frequency 
deviation,  h  =  2f^T,  where  T  is  the  bit  time.  The  signal  bandwidth  is  propor¬ 
tional  to  h.  At  the  receiver,  an  I.F.  (intermediate  frequency)  filter  excludes 
unwanted  receiver  conversion  products  and  controls  the  amount  of  noise  admitted. 
If  the  bandwidth,  Wjp,  of  this  filter  is  large,  the  received  FM  waveform 
passes  through  to  the  detector  undistorted.  If  the  bandwidth  is  decreased, 
less  noise  is  admitted  but  the  resulting  distortion  of  the  waveform  gives 
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rise  to  intersymbol  interference.  Thus  there  are  optimum  values  of  h  and 
the  product  WjpT  which  are  traded  off  against  one  another.  Without  hopping 
or  jamming,  it  has  been  shown  the  h  =  0.7  and  WjpT  =  1.0  are  best  values  [3]. 
With  worst-case  partial-band  jamming  but  neglecting  thermal  noise,  it  has 
been  asserted  that  h  =  0.6  and  WjpT  =  0.75  are  best  values  [22]. 

1.2  RATIONALE  FOR  THE  STUDY 

We  have  mentioned  previously  that  detailed  analyses  of  hopped  CPFSK 
systems  under  various  kinds  of  jamming  and  also  system  or  thermal  noise  are 
not  available  in  the  current  literature.  Because  of  the  potential  advantages 
offered  by  CPFSK  over  MFSK  and  because  many  current  hopped  and  unhopped  radio 
systems  employ  some  form  of  CPFSK,  it  is  important  to  expand  the  knowledge  of 
the  performance  that  can  be  expected  from  FH/CPFSK  systems,  and  what  parameter 
values  are  optimum. 

Noncoherent  detection  of  CPFSK  can  be  performed  using  a  limiter-discrimi¬ 
nator  with  integrate-and-dump  filtering,  or  by  differential  detection.  In 
either  case,  a  difference  L‘  in  the  total  I.F.  waveform  phase  is  extracted 
which,  in  the  absence  of  noise  and/or  distortion,  would  convey  the  transmitted 
binary  information.  Normally,  that  is,  without  repetition  or  diversity,  a 
hard  decision  is  made  based  on  the  sign  of  A4.  Typically  A4  is  corrupted  by 
noise  projected  onto  the  phase  of  the  signal  by  the  action  of  the  bandpass 
limiter  which  precedes  the  discriminator.  Now,  since  the  phase  of  a  carrier 


is  inherently  ambiguous  (modulo  2tt),  it  can  be  written 
A 4  =  (A4 )  mod2^  +  2^N, 


(1.2-1) 


where  N  is  the  net  number  of  positive  2r  phase  rotations  or  "clicks"  in  the 
bit  interval;  usually  N  is  taken  to  be  a  negative  random  (Poisson)  integer 
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if  the  data  gives  At  >  0,  and  positive  if  A4>  <  0.  For  high  carrier-to-noise 
ratios  (CNR),  the  probability  that  N  is  nonzero  becomes  very  small.  It  has 
been  argued  then  that  the  limiter-discriminator  detection,  in  effect,  limits 
or  clips  the  phase  noise,  and  therefore  that  jamming  effects  in  an  FH/CPFSK 
system  will  be  suppressed.  Our  study  tests  this  hypothesis  by  postulating 
and  evaluating  an  FH/CPFSK  system  with  soft-decision  diversity  combining 
of  received  ai>  samples,  implementing  the  decision  rule 

bit  =  1 

SUM  1  ( A ) .  }  *  0.  (1.2-2) 

k  k  < 

K  bit  =  0 

In  our  analysis  and  computations  forming  the  body  of  this  report,  we 
follow  a  rigorous  treatment  of  the  system,  including  background  noiie, 
intersymbol  interference,  and  FM  noise  clicks.  We  consider  linear  combining 
of  receiver  samples  for  both  1 imi ter-discriminator  and  differential  detection 
types  of  CPFSK  receivers,  since  both  are  used  in  current  tactical  FH/CPFSK 
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1.3  SYSTEM  STUDIED 

In  view  of  the  importance  of  the  FH/CPFSK  waveform  to  tactical  military 
communications,  and  of  the  potential  for  diversity  improvements  in  hopping 
system  performance  against  jamming,  we  have  undertaken  the  studies  summarized 
in  this  report.  The  following  material  briefly  describes  the  system  studied. 
1.3.1  Transmission  Scheme 

Our  studies  concern  the  jammed  performance  of  hopped  binary  FM  communica¬ 
tions,  particularly  under  the  assumption  of  partial-band  noise  jamming  and  the 
use  of  (time)  diversity  to  mitigate  the  jamming.  Figure  1.3-1  gives  a  block 
diagram  of  the  transmission  scheme  for  the  system.  Binary  data,  after  error- 
control  coding,  is  to  be  transmitted  using  slow-frequency-hopping  digital  FM, 
or  CPFSK.  The  coded  symbols  are  to  be  repeated  on  L  different  hops  in  order 
to  increase  the  likelihood  that  some  of  the  symbols  are  free  of  jamming.  The 
figure  suggests  one  of  many  possible  ways  to  accomplish  this  objective. 

According  to  the  version  shown  in  the  figure,  the  coded  symbols  are  first  read 
into  a  Q-bit  shift  register  (Q-symbol  buffer),  where  Q  is  the  number  of  symbols 
that  can  be  transmitted  in  one  hop  period.  For  example,  if  the  channels  allotted 
to  the  system  support  20  kbps  digital  FM  signalling,  and  the  hop  rate  is 
100  hops/sec,  then  Q  could  be  200. 

When  the  Q  symbols  have  all  been  generated  at  rate  Rc  and  stored  in  the 
input  buffer,  they  are  then  transferred  to  a  second  (output)  buffer.  The 
transmitter  logic  reads  this  buffer  L  times  at  the  rate  LRC,  and  this  stream 
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of  data  "chips"  is  used  to  frequency-modulate  the  selected  carrier 
frequency,  which  is  changed  (hopped)  to  a  new,  pseudorandomly-selected  value 
after  Q  chips  have  been  transmitted.  In  this  manner,  L  copies  of  the  Q-symbol 
sequence  have  been  transmitted  on  L  different,  successive  hops.  Although 
Figure  1.3-1  suggests  that  the  Q  symbols  are  repeated  in  the  same  order  on 
each  hop,  it  is  of  course  possible  with  a  more  sophisticated  system  to  permute 
the  symbols  or  otherwise  scramble  them  so  that  the  ordering  of  the  symbols  is 
different  on  each  hop. 

We  note  that  certain  fundamental  relationships  exist  among  the  digital 
rates  at  various  points  in  the  transmission  logic,  and  among  the  energies  in  the 
transmitted  waveform  which  correspond  to  each  rate.  The  symbols  actually  trans¬ 
mitted  are  keyed  at  the  rate  Rrf  which,  as  we  have  already  noted,  is  a  basic 
specification  of  the  communications  channels  being  used  for  hopping.  The 
original  symbol  rate  is  Rc  =  Rd/L,  on  account  of  the  repetitions.  Viewed  another 
way,  the  energy  transmitted  per  chip  is  the  fraction  1/L  of  the  coded  symbol 
energy.  If  the  error  control  code  rate  is  r,  then  the  original  bit  rate  is 

Rb  =  Rcr  =  r  Rd/L’  (1-3-1) 

For  example,  for  =  20  kbps,  r  =  1/2,  and  L=5,  Rc  =  4  kbps;  between  the 
coding  and  the  repetitions,  the  bit  energy  gets  split  up  into  10  pieces  in  this 
example. 

1.3.2  Reception  scheme 

It  is  obvious  that  synchronization  and  timing  are  especially  important 
to  the  type  of  waveform  we  are  considering.  The  receiver,  diagrammed  in 
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Figure  1.3-2,  must  accurately  synthesize  local  oscillator  frequencies  in  order 
to  dehop  the  signal  at  the  proper  times,  and  must  then  be  capable  of  sampling 
the  demodulator  output  at  the  ends  of  each  chip  interval.  These  samples,  QL 
of  them  for  L  complete  hops,  are  buffered  so  that  the  receiver  logic  can  in  some 
manner  combine  the  L  chips  belonging  to  a  particular  code  symbol. 

As  the  figure  suggests,  binary  decisions  can  be  made  on  each  chip  as  it 
is  received.  The  resulting  logic  and  buffering  for  such  a  "hard  decision" 
procedure  is  simpler  than  that  required  for  a  "soft  decision"  scheme,  which 
involves  A/D  conversion  of  the  samples  and  storing  the  resulting  QL  multi-bit 
words,  one  per  chip.  With  either  hard  or  soft  chip  processing,  after  diversity 
combining  there  is  the  option  to  do  binary  decoding  of  code  symbol  decisions,  or 
soft-decision  decoding  of  code  symbol  metrics  produced  by  the  combining. 
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1.4  SUMMARY  AND  RECOMMENDATIONS 
1.4.1  Main  Conclusion 


Our  work  thus  far  on  FH/CPFSK  systems  with  diversity,  documented  in 
this  report,  has  determined  that  linear  combining  or  summing  of  the  diversity 
components  (chips)  of  each  bit  does  not  yield  an  improvement  in  uncoded 
system  performance  against  worst-case  partial-band  noise  jamming.  This 
conclusion  holds  for  both  limiter-discriminator  and  differential  detection 
techniques . 

1.4.2  Achievements 

In  developing  analytical  and  computational  methods  for  the  calculations 
which  produced  this  conclusion,  we  have  made  considerable  advances  in  the 
methodology  for  evaluating  FH/CPFSK  and  CPFSK  systems.  These  include: 

(a)  Independent  derivations  of  both  integral  and  series  forms  of  the 
differential  phase  probability  density  function,  showing  explicitly  the 
effect  of  taking  the  phase  to  be  modulo  2r  . 

(b)  Explicit  derivations  and  formulas  for  intersymbol  interference-related 
SNR  and  differential  phase  parameters  for  all  eight  of  the  possible  adjacent 
bit  data  patterns,  and  arbitrary  I.  F.  filter  transfer  functions. 

(c)  An  independent  and  detailed  derivation  of  FM  noise  click  rates  and 
average  number  of  positive  clicks,  and  development  of  the  notion  of  "signifi¬ 
cant  clicks"  in  computation  of  the  error  probability. 

(d)  A  general  formulation  of  the  L  hop/bit  FH/CPFSK  jammed  error  probability 
valid  for  any  diversity  combining  technique. 

(e)  An  original  derivation  of  the  diversity  sum  jammed  error  probability 
for  FH/CPFSK  using  differential  detection. 

(f)  Comparison  of  exact  and  simplified  performance  calculations  (the 
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exact  require  much  computer  time  and  effort;  the  simplified  can  be  done  on 
a  programmable  calculator). 

(g)  Development  of  a  computer  program  based  on  the  discrete  Fourier 
transform  for  evaluating  the  probability  of  error,  which  is  useful  for  all 
sum  metrics. 

(h)  Demonstration  that,  if  perfect  side  information  is  used,  there  is 
a  diversity  gain  for  FH/CPFSK. 

1.4.3  Recommendations 

The  diversity  sum  method  for  combining  the  L  chips  for  a  given  FH/CPFSK 
bit  can  be  regarded  as  implementing  a  form  of  soft-decision  metric.  Note 
that  the  summing  operation  does  not  utilize  any  side  information  on  whether 
a  particular  chip  is  jammed,  or  how  strongly  it  is  jammed.  Since  there  is 
no  mechanism  for  excluding  or  otherwise  treating  jammed  chips  differently, 
when  one  or  more  chips  are  subject  to  jamming,  the  entire  metric  is  corrupted, 
and  this  accounts  for  the  failure  of  the  linear  combining  (simple  sum)  metric 
to  provide  a  diversity  improvement. 

Now,  it  is  a  fact  that  if  instead  of  soft-decision  combining  of  the 
chips,  we  combine  hard  decisions  on  each  chip,  for  high  E^/Nq  (15  dB  or  more) 
a  diversity  gain  results,  in  the  sense  that  for  a  particular  J/S  ratio  L  >  1 
may  yield  a  lower  bit  error  probability.  The  improvement  is  due  to  the  hard 
decisions'  limiting  a  jammed  chip  to  "one  vote"  in  the  sum  of  hard  decisions. 
Since  no  side  information  is  required,  this  simple  ECCM  scheme  is  very 
attractive,  except  for  the  fact  that  noncoherent  combining  losses  are  high 
for  hard-decision  metrics,  relative  to  soft  decisions. 

It  can  be  shown  that  a  "perfect  side  information"  soft-decision  scheme 
which  includes  only  unjammed  chips  in  the  sum  can  provide  greater  diversity 
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gain  than  the  hard-decision  combining  strategy.  Therefore,  some  practical 
method  which  combines  analog  chip  samples  or  soft  decisions  is  likely  to 
exist  for  improving  FH/CPFSK  performance  in  partial-band  jamming  more  than 
the  hard-decision  metric  does.  Such  metrics  have  been  found  for  noncoherent 


FH/BFSK  [23,24],  including  a  "self-normalizing"  technique  which  does  not 
require  side  information.  Therefore,  we  recommend  that  further  studies  of 
the  type  given  in  this  report  be  conducted  of  soft-decision  diversity 
combining  metrics  for  FH/CPFSK  which  have  the  potential  for  improved  perform¬ 
ance  in  jamming. 
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2.0  MODELLING  AND  PRELIMINARY  ANALYSES 

In  this  section,  we  describe  the  modelling  assumptions  used  to  analyze 
the  performance  of  frequency-hopped  continuous  phase  frequency-shift  keying 
(FH/CPFSK)  in  the  presence  of  thermal  noise  and  partial-band  noise  jamming. 
Limiter-discriminator  (LD)  detection  is  assumed.  Certain  preliminary  analyses 
are  conducted  to  predict  the  effects  of  intersymbol  interference  (ISI)  on  the 
detected  waveform,  and  to  derive  the  distribution  of  the  noise  projected  onto 
the  signal  phase. 

2.1  DESCRIPTION  OF  SYSTEM  WAVEFORMS 

For  discussion  purposes,  we  consider  first  the  simplified  system  model 
depicted  in  Figure  2.1-1. 

2.1.1  Overview  of  System  Operation 

A  binary  data  source  outputting  the  sequence  { dk)  modulates  an  FM 
transmitter  with  the  NRZ  bipolar  data  waveform  d(t),  where 


d(t)  =  ^  dk  p(t-kT) ,  dk  =  +  1, 


(2.1-1) 


and  the  pulse  function  p(t)  is  assumed  to  be  rectangular: 


p(t)  =  u( t)-u(t-T) 


=  rect(t-T/2) , 


with  u(t)  the  unit  step  function.  The  interval  T  is  the  bit  duration,  so 
that  the  data  rate  is  Rfa  =  1/T.  Although  much  attention  is  being  directed 


(2.1-2) 


KfivHKC 


» 


I 

: 

i 


currently  at  reducing  the  transmitted  bandwidth  of  binary  FM  signals  by 
employing  non-rectangular  p(t)  and  by  correlative  coding  of  the  data,  here 
we  will  assume  the  p(t)  as  given  and  will  treat  the  d^  bit  values  as  having 
been  generated  independently.  (A  more  detailed  discussion  of  the  generation 
of  the  data  will  be  presented  below.) 

The  resulting  modulated  carrier  frequency  is  given  by 


f(t)  =  fQ  +  fd  •  d(t). 


and  the  commonly  accepted  measure  of  relative  frequency  deviation  is  the 
modualtion  index 


h  =  2fd/Rb  =  2fdT-  <2- 

Thus  the  emitted  waveform  (neglecting  any  post-modulator  filters  or  distortion, 
and  without  frequency  hopping  of  the  carrier  fg)  is 


s(t)  =  const  >•  cos  2r fpt  +  2rfd 


di  d(t ) 


Ideally,  the  frequency  hopping  and  dehopping  operations  shown  in  Figure 
2.1-1  are  in  perfect  synchronism  (with  the  propagation  delay  accounted  for),  and 
are  therefore  "transparent"  in  the  sense  that  whether  the  carrier  is  fixed  or 


hopped  should  not  affect  system  operation. 
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Against  a  noise  jammer  occupying  a  portion  of  the  RF  band  over  which  the 
system  is  hopping,  the  hopping  and  dehopping  result  in  the  following  total 
waveform  prior  to  intermediate  frequency  (I.F.)  filtering,  overlooking  harmonics 
to  be  rejected  by  that  filtering: 


r(t)  =  A  cos{2irf0t  +  eQ  +  em( t) }  +  n(t). 


(2.1-6; 


In  this  formulation,  we  show  an  undistorted  signal  term;  ~q  is  a  random  phase 
offset  left  over  from  the  downconversion  and  dehopping,  and  for  convenience 
we  define 


-m(t)  s  2rf . 
m  d 


The  noise  term  n(t)  is  given  by 


n(t)  = 


d-r  d(  x) . 


n0(t),  hop  was  not  jammed 


nQ(t)  +  nj(t),  hop  was  jammed. 


(2.1-7) 


(2.1-8) 


The  background  and  receiver  additive  bandpass  noise  term  n^(t)  and  the 
jamming  bandpass  noise  term  n^(t)  are  assumed  to  have  flat  spectra  going  into 
the  I.F.  filter,  with  two-sided  spectrum  levels  Nq/2  and  Nqj/2,  respectively. 

Conceptually,  the  received  signal  -  with  more  or  less  noise,  depending 
on  whether  the  signal  was  jammed  or  not  on  a  particular  hop  -  is  filtered  and 
demodulated  to  recover  an  estimate  d(t)  of  the  data  waveform  d(t)  sent  originally 
by  the  transmitter.  The  information  in  d(t)  is  then  extracted  by  receiver 

processing. 
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Now  we  must  go  into  considerably  more  detail  to  describe  the  structure 
of  the  data  sequence  and  its  relation  to  the  hopping  scheme.  Following  that, 
we  discuss  the  distortion  effects  of  receiver  filtering  and  describe  the 
receiver  processing  in  further  detail. 

2.1.2  The  I . F ,  Filter  Output,  Signal  Term 

The  I.F.  filtering  operation  is  intended  to  reject  all  but  the  selected 
signal.  Its  finite  bandwidth  is  taken  to  be  its  noise  bandwidth  WJf:  =  2Bg, 
where  the  bandpass  filter  transfer  function  H(f)  is  modelled  in  terms  of  a 
lowpass  filter  transfer  function  Hg(f)  by  the  relation 

H(f)  =  Vf“f0}  +  H0(f+f0);  (2.1-9a) 

that  is,  its  impulse  response  is  considered  to  be 

h(t)  =  2hg( t)  cos  .» q t ,  ( 2 . 1  -9 b ) 

and  hg(t)  is  the  lowpass  filter  impulse  response  corresponding  to  Hg(f). 

Accordingly,  the  noise  bandwidth  Bg  is  given  by 


df  H0(f)  ‘  /  |Hg(0), 


(2.1-10) 


For  example,  if  the  lowpass  filter  has  the  Gaussian  shape 


H0(f)  =  e 


-rf^/8B ; 


(2.1-11) 


then  the  3  dB  bandwidth  Bg  is  0.9294  times  the  noise  bandwidth  Bg.  For.  an 
n-pole  Butterworth  filter  with  transfer  function 
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A  J  dt  h0(t-T)  cosem(i)  {  cos(w0t+e0)+  cos[w0(2T-t  )+?Q] } 


-A J  dr  h^t-r)  sinem(T){sin(wot+©Q)+sin[a)o(2T-t)+eQ]} 

—  oo 

=  A  cos (oj0 t+e^) J  dT  hg(t-T)  cosem(-r) 

-  A  sin(u)ot+6Q)  J  di  hgU-i)  sine^-r) 


=  a(t)  A  cos[ui0  t+e0+®(t)],  fn  »  B, 


(2.1-13) 


where 


a- (t)  =  [Mt)  *  cosi_(t)]2  +  [ hn( t)  *  sine  (t)]2 


(2.1-14a) 


(t)  =  tan” ~ [nQ( t)  *  sinem(t)]/[h0(t)  *  cosem(t)]}  .  (2.1-14b) 


From  this  development  we  see  that  the  filter  distorts  the  signal  phase  waveform 


from  -  (t)  to  : ( t) ,  and  induces  amplitude  modulation  a(t). 


For  most  cases  of  practical  interest,  it  is  sufficient  to  consider  inter¬ 


symbol  interference  effects  (i.e.,  the  overlapping  of  filter  responses  from 
different  bit  intervals)  due  to  immediately  adjacent  bits  [1].  Therefore,  in 


what  follows  we  consider  bit  patterns  which  are  the  periodic  extensions  of  the 


patterns 


111,  000  (all  one's  or  zeros); 

010,  101  (alternating  one's  and  zeros); 


(2.1-1 5a ) 


(2.1- 15b ) 


0110,  1001,  1100,  0011. 


(2.1-1 5c ) 


The  "present  bit"  in  these  sequences  is  indicated  by  the  underlining.  The  patterns 
in  (2.1-15)  were  chosen  because  they  generate  the  eight  possible  3-bit  patterns 


in  a  very  simple  manner  and  can  be  analyzed  easily. 
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Using  the  steady-state  filtering  approach  of  Tjhung  and  Wittke  [2]  and 
Pawula  [3],  we  recognize  that  if  the  Fourier  series  for  the  periodic  extension 
of  the  i : th  patterns  yields  (assuming  evenness  of  about  t  =  0) 


co 

sin^U)  =  Y.  a^k  cos(2*kfpt) 


(2.1 - 16a ) 


( 1 )  r—  ( } 

COS5m  =  L  6  k  C0S(2*kfptJ 


(2. l-16b) 


where  f  is  the  appropriate  fundamental  frequency  of  8^' ' ( t) ,  then  the  responses 
of  the  lowpass  filter  hg(t)  to  these  components  are 


u^t)  =  hQ(t)*sine^(t) 


oo 

=  Z  |H0(kfp)lakl}  c°sC2Trkf  t  -  B(kf  )] 
k=l  h  h 


( 2 . 1- 17a ) 


and  vi(t)  =  hQ(t)  *  cose^(t) 


uu 

=  X  lHo(kV  1(4°  cosl>kf  t  -  B(kf  )] 


(2.1-17b) 


where  B(f)  is  the  filter  phase  delay.  For  convenience,  we  shall  employ  the 
alternative  notations 

u(t;pattern  i)  h  u.(t) 
v(t;pattern  i)  s  v^t)  . 


(2.1-17C) 
(2. l-17d) 


Note  from  (2.1-15)  that 


tani(t)  =  u(t)/v(t) 


(2.1-18) 


for  each  pattern, 
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For  patterns  111  and  000,  the  signal  is  a  pure  sinusoid,  so  that 
(neglecting  any  filter  delay)  the  quadrature  components  of  the  signal  are 

u(t;lU)  =~u(t;000)  =  aQ  sin(irht/T)  (2.1-19a) 

and 

v(t;Ul)  =  v(t;000)  =  aQ  cos(nht/T).  ( 2 . 1-  19b ) 

That  is,  for  these  two  patterns,  a(t)  =  and  4> ( t )  =  ±ht/T,  where 

a0  =  lH0(fd>l  •  (2.1-190 

For  patterns  OK)  and  101,  the  original  frequency  modulation  is  a  ±f^ 
squarewave  with  period  2T,  so  that  *m(t)  is  a  bipolar  triangular  wave  with 
amplitude  Th/2  and  period  2T.  For  pattern  010,  the  triangular  wave's  positive 
peak  occurs  at  t  =  0,  using  the  convention  that  d(t)  =  dk  for  (k-l)T  <  t  <  kT. 
Expanding  sin  e  ( t)  and  cos  em(t)  in  Fourier  series  gives 


sine  (t;010)  =  -sineJtUOl) 


cos(y) 


V  cosr(2k-l)Tt/Tl 
(2k-l)2-h2 


(2.1~20a) 


and 


cos?m(t;0_10)  =  cos9m(t;l£l) 


2_ 

*h 


sin(^y) 


l 


k=l 


cos(2kTt/T) 

( 2k )2- h2  . 


(2.1 -20b ) 


Assuming  that  the  I.F.  filter  passes  only  harmonics  of  these  Fourier  series 
up  to  f  =  1/T,  we  find  that 


'M 


u  ( t ;  OK) )  =  -u(t;101) 


it  C0S(1I2  ^!H0^2T^I  7T2  cos^T 

1-h 


B<2T» 


(2.1-21a) 


v(t;010)  =  v{  t ;  10.1 ) 

7h  sin(y) 


1  ‘  fV  lHo{T)|cos[1r  ‘  b(t)]  » 


(2.1-21b) 


where  the  phase  delay  B(f)  is  taken  to  be  zero  if  a  Gaussian-shaped  filter  is 
assumed.  Figure  2.1-2  illustrates  the  various  waveforms  associated  with  this 
type  of  pattern. 

For  bit  patterns  0110  and  1001*  the  original  frequency  modulation  is  a 

squarewave  with  period  4T,  so  that  em(t)  is  a  bipolar  triangular  wave  with 

amplitude  irh  and  period  4T.  For  pattern  0110,  the  triangular  wave's  positive 

peak  occurs  at  t  =  T.  By  analogy  with  (2.1-20),  the  Fourier  series 

expansions  for  sine  (t)  and  cose  (t)  are 

m  m 


sine  (t;H00)  =  -sine  (t;0011) 


Si  cos(„h)  y  cos[(2lc-l)i.(t- 

"  k4j  (2k-l)2-4h2 


(2.1-22,) 


and  cosem(  t;lj_00)  =  cosem(  t ; 001 1 ) 


sinUl 

irh 


1  -  2h‘ 


cos[  kiT  ( t-T  V/T] 


(2.1 - 22b ) 


Assuming  the  filter  rejects  harmonics  with  f  >  1/T,  we  find  that  the  signal 
quadrature  components  for  these  patterns  are 
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u(t;U00)  =  -u(t;0011) 


—  cos(Tih)  -1  lHn(Tf)|sin[n/2T-B(l/4T)] 
71  1  -4ir  u 


■K  I  sin[3*t/2T-B(3/4T)] 


(2.1-23a) 


v( t ; 1 100 )  =  v(t;0011) 


1  +^|H0(^f)|cos[irt/T-B(l/2T)] 


~2  I  Hq ( y)  I  COS[2nt/T-B(l/T)] 


(2 . 1 - 23b ) 


Figure  2.1-3  illustrates  the  various  waveforms  associated  with  this  type  of  pattern. 

Recognizing  that  the  time  extensions  of  the  patterns  0110  and  1001  are  time 
shifted  versions  of  those  for  1100  and  0011,  we  can  immediately  write 


u(t;0U0)  =  -u(t;1001) 

=  u(t+T;lJ_00) 

v(t;0U0)  =  v  ( t ;  100 1 ) 

=  v(t+T;1100) 


(2,l-24a) 


(2. l-24b) 


A  summary  of  these  components  which  determine  4>(t)  is  given  in  Table 
2.1-1,  assuming  no  filter  delay,  with  example  values  as  listed  in  Table  2.1-2 
for  h  =  0.7  and  D  =  WjpT  =  1.0. 


n 

$ 

§ 

ffi 

?■< 

Si- 

•J 

MC 

r$l 

S 


.*  «  * 

» 

SL 


vv,v;/v.y;.- v; 


TABLE  2.1-2 


EXAMPLE  VALUES  FOR  SIGNAL  PARAMETERS  DESCRIBED  IN  TABLE  2.1-1 
Case:  h  =  0.7,  D  =  WjpT  =  1.0,  Gaussian  filter 


Parameter 


Value 

0.82496 


0.67523 


0.20788 

0.90649 

0.41330 


0.79339  a1  =  0.53572 
0.81033 

0.22625  a,  =  0.047032 


1.0914  a3  =  0.98935 
-0.14883  a4  =  -0.061511 
0.36788 

0.70691  al  =  0.47732 
0.10271  a2  =  0.021352 


WWW 


ill 
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2.1.3  Limiter  Output  and  Total  Phase 

A  simplified  receiver  block  diagram  was  given  previously  as  Figure  2.1-1. 
As  shown  in  Figure  2.1-4  the  digital  FM  demodulator  consists  of  a  limiter 
followed  by  a  discriminator  and  an  integrate-and-dump  filter,  whose  output 
is  sampled  to  yield  the  demodulated  data  sequence. 

The  purpose  of  the  limiter,  assumed  to  be  an  ideal  bandpass  limiter,  is 
to  remove  any  amplitude  modulation  on  the  I.F.  signal.  Its  output,  y(t),  is  a 
constant  amplitude  sinusoid  with  the  same  total  phase  (including  modulation) 
as  the  I.F.  waveform: 


y(t)  =  constant  ■  cos:~;t  +  Phase  [x( t ) ] ; 


The  total  phase  ;(t)  of  x(t)  is  found  from  the  following  development: 


(2.1-25; 


x(t)  =  [a(t)A  cos : ( t)  +  nc(t)]  cos(~ct+->0) 

-  [a(t)A  s i n : ( t )  +  n s ( t ) ]  sin(^ot+cQ) 


=  Env[x( t)]  •  cos[-;t  +  -n  +  c(t)  +  n(t)]. 


where  r(t)  is  a  phase  noise  term,  and 


tan[ :( t)  +  n(  t)]  -  tan?(t) 


a(t)A  sin-(t)  +  ns(t)  u(t)  +  n$(t) 
a(t)A  cos:(t)  +  n(t)  '  v(t)  +  n.(t) 


This  development  uses  the  Rician  decomposition  of  the  bandpass  noise, 
referenced  to  the  carrier  frequency  and  phase: 


(2.1-26a; 


(2 . 1 -26b ) 


raws 


r . . 

j  J.  S.  LEE  AS 

SOCIATES,  INC. 

[ 

n(t)  =  nc(t)  cos(.**0 1+?0)  -  n$(t)  sin(u>0  t+eQ) . 

Recall  that  this  noise  may  include  jamming  noise.  At  any  instant,  the 
quadrature  components  nc(t)  and  ng(t)  are  independent  zero-mean  Gaussian 
random  variables  with  equal  variances 

NgWjp,  hop  not  jammed 

= 

( No+NOJ )WI F *  h0p  jammed- 
Note  that  we  can  write  the  phase  noise  term  as 


(2.1-26c) 


(2.1-27) 


r,  (t)  =  tan 


_1  r  v$(t)cos;(t)  -  vc(t)sin$(t) 

. v'2o(tT"  +  v  (t)sin<^(t)  +  v  (t)coso(t) 

5  v 


(2.1-28a) 


in  which  v  (t)  and  v  (t)  are 

o  v*. 


v$(t)  -  n$(t)/c  ,  vc(t)  =  nc(t)/o  , 


(2.1- 28b ) 


unit-variance  Gaussian  random  variables,  and 

s(t)  *  a2(t)A2/2c2  =  [u2(t)+v2(t)]/2o2  (2.1-28c) 

is  the  SNR  (carrier-to-noise  ratio)  with  time  variation  due  to  the  I.F. 
filter-induced  amplitude  modulation.  A  further  simplification  results  from 
using  a  rotational  transformation  of  noise  variables  to  write 

1  r  1 

_  ■*.,_  •*•  •>  .  lo  i  on-.1! 


r(t)  =  tan’ 

we  recognize  that  since 


£p(t)  +  v'r(t)J 


(2.1- 29a ) 


'w'*  . 

v"[vyB> 


S’1 

5 

&i 


■jO\  >: 


•• .  v- 
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i  *  \ ^  »  -  %2 

(vs)  +  (vc) 

/  .  2  2 
=  (v$  cos  $  -  vc  sin<j>)  +  (\>s  s1n«^  +v£  cos<f>) 

=  v 2  +  V2  , 

S  C  ’ 

the  distribution  of  the  rotated  noise  components  v  and  v  is  identical  to 

c  c 

that  of  the  actual  components  vc  and  vg.  Thus  at  a  given  time  instant,  the 
phase  noise  termn(t)  is  independent  of  the  signal  phase  and  additive  to  it. 
It  is  easily  shown  that  the  probability  density  function  (pdf)  for  n(t)  is 
(see  [4],  chapter  9,  e.g.) 


(2.1-29b) 


I 

U'lj 

ra 

I 


pn(a)’fr(  dx  x  e^-  **  -•»« 
J  n 


(2.1-30a)  $ 


9 


*  ^  e"PSin  “  cosa  Q(-»/?P  cosa),  |a|<*  • 


» 


and  therefore  the  total  phase  pdf  is 


(2. i-30b) 

•.v 

■-V 


p.(a)  =  p  (a -♦) 


=  +  e"psin^a*^  cos  (a  —  4> )  Q[-/2P  cos  (a  -<t>)  ] ,  |  o-4 1  <n .  (2.1-30c) 


In  these  expressions,  Q( •)  is  defined  as 


=  i  erfc  (x//2). 


BS 


(2. l-30d)  K3 
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2.1.4  Discriminator  and  Baseband  Outputs 

The  discriminator  output  extracts  2*  times  the  instantaneous  frequency 
deviation  from  the  carrier  fg.  This  quantity  is 

z(t)  =  2*[f(t)-f0]  =  4*(t)  +  n(t)  =  *(t).  (2.  l-31a) 

The  baseband  filter,  assumed  to  be  of  the  integrate-and-dump  type,  operates  on 
z(t)  to  produce  the  differential  phase  at  sample  time  t$: 


Although  we  have  written  the  total  differential  phase  A4>  as  the  sum  of  a 
signal  differential  phase  term  and  a  noise  differential  phase  term,  in  general 
the  "differential  phase  noise"  term  An  is  not  additive  and  independent  of  a$, 
but  rather  depends  upon  A$..  Further  discussion  of  phase  noise,  including 
"FM  clicks",  is  given  in  Sections  2.1.5  and  2.1.6. 

Without  noise,  the  differential  phase  output  of  the  digital  FM  demodulator 
is 


-1  U<V  -1 

"  tan  [vrg'J  "  tan  v(ts-T)  *  (2. l-32a) 

where  the  principal  values  of  the  arctangents  can  be  used  if  the  trajectory 
of  *(t)  =  tan-1  [u/v]  is  such  that  |a$|  <  tr  .  Otherwise,  there  is  an  inherent 


2:r-radian  ambiguity  in  the  arctangent.  For  the  signal  only,  these  requirements 
are  satisfied,  and  the  resolution  of  arctangent  ambiguities  can  be  successfully 


"u>- 
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accomplished  by  rewriting  (2.1-31a),  using  the  relation 


9i  -  s2  =  tan-1  {  tan(e1-e2 )} 


tan  j  — — — — — v  , 

(  1  +  tane1tane2  J 

and  identifying  6j  and  s2  to  be  the  first  and  second  terms,  respectively  of 
(2.1-32a),  it  is  straightforward  to  find  that 


_ j  j  tan  e j  ■  tan  02 


Ad  =  tan 


-i  f  u(Vv<VT)  -  u(VT)v(ts> 


ults)u(ts- 


v(ts)v(ts- 


(2. 1 - 32b ) 


A  plot  of  Ad(t)  for  the  various  data  sequences,  commonly  called  an  "eye 
pattern",  is  shown  in  Figure  2.1-5  for  h  =  0.7  and  D  =  1.0.  Note  the  "closing" 
of  the  "eye"  for  the  alternating  bit  sequences,  due  to  the  filtering  in  the 
receiver,  and  the  dependence  of  the  Ad  values  at  the  bit  times  upon  the 
sequence. 

The  value  of  A?(t)  at  t  =  0  represents  the  sampled  data  value  as  recovered 
by  the  demodulator.  From  Table  2.1-1  we  can  calculate  the  sample  amplitudes 
observed  in  Figure  2.1-5  in  terms  of  the  patterns  and  filter  parameters.  These 
are  summarized  in  Table  2.1-3.  From  the  figure  and  table  we  observe  the  effect 
that  the  I.F.  filter  has  on  the  data  output.  Ideally,  the  differential  phase 
is  +irh,  as  it  is  for  the  all-one's  or  all-zeros  patterns.  For  the  other  patterns, 
the  phase  distortion  causes  intersymbol  interference,  as  the  values  of  adjacent 
bits  affect  that  of  the  current  one,  with  the  worst-case  being  alternating 
one's  and  zeros.  However,  as  Table  2.1-3  demonstrates,  the  average  of  the 
time-varying  SNR  is  higher  for  the  alternating  patterns,  since  the  instantaneous 
frequency  for  these  patterns  is,  on  the  average,  nearer  to  the  filter  center 
frequency  than  is  that  for  the  non-alternating  patterns. 
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2.1.5  Differential  Dhase  Hi str’bution 


As  we  have  seen,  the  digital  FM  receiver  output  produces  the  differential 


phase  M(t),  where  the  difference  is  between  the  two  values  of  the  signal -plus- 


bandpass  noise  total  phase  occurring  at  the  sampling  time  and  one  bit  period 


earlier.  The  probability  distribution  of  the  differential  phase  is  found 


by  starting  with  the  joint  distribution  of  the  additive  noise  components.  Let 


the  noise  quadrature  components  be  given  by 


ncU)  =  n!  nc(t-T)  e  n3 


n$(t)  =  n2  n$(t-T)  e  n4 


( 2 . 1 -33a ) 


If  the  autocorrelation  function  for  the  bandpass  noise  is  given  by 


Rn(0  =  ~  [r( t)  cos~ct  -  >(t)  sin.0x]. 


(2. 1 - 33b ) 


then  the  column  vector  ji  =  (n^  n2>  n3>  n^)  is  a  zero-mean  multivariate 


Gaussian  random  vector  with  covariance  matrix 


r  x 


=  ozC, 


(2.  l-33c' 


where 


r  e  r(T)  and  >.  e  \  (T) . 
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For  a  flat  noise  spectrum  going  into  the  receiver  I.F.  filter,  the  auto¬ 
correlation  function  Rn(t)  is  determined  by  the  shape  of  the  filter.  For 
example,  if  the  filter  frequency  characteristic  (passband)  is  symmetric  about 
the  signal  center  frequency  fg,  then  a(t)  is  identically  zero.  Using  the 
Gaussian-shaped  filter  introduced  previously  in  equation  (2.1-11),  the  noise 
spectrum  is 


“r  ( f“f0 )  2/WIF  -1T(f+V2/WiF 

e  +  e 


(2.1-34a) 


for  which  the  autocorrelation  function  is 


-r(W  t)^ 

Rn(')  *  N0Wjf  e  IF  cos-'qt , 


( 2 . 1- 34b ) 


so  that  for  this  filter 
■*2 


r  -  e 


-~D 


=  0, 


=  N0WIF 


(2.1-34c) 


in  (2.1 - 33b ) ,  using  D  =  WjpT. 

The  multivariate  probability  density  function  for  the  noise  vector _n  is 


= 


ir-^bet  I  ) 


-1 


exp 


1  T 


-~2  o_  E 


(2.1-35a) 


If  we  define  the  signal-plus-noise  vector  as  x,  with  signal  component  Xg  given 
by 

r 

Aaj  cosh 


^0 


Aa  j  s  i  n  i 
Aa2  cosh 


Aa  si  ncr  2 

2 


(2.1 - 35b ) 
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then  the  pdf  for  _x  is 


Px(iL)  =  Pn(£-x<)); 


for  convenience  we  can  normalize  each  noise  component  by  a  to  get  unit-variance 
random  variables,  with  the  normalized  signal  component  vector  (mean)  then 
expressible  as 


1  =  v~ 


/2  COS:  j 

/FT]  S  i  n  c  i 
/FT]  cosc2 
\_/F72  sin ?,2 J 


(2. l-35d) 


with  :■  and  being  the  SNR  at  times  t  and  t-T,  respectively.  The  normalized 
pdf  then  is 


-1 


PXU)  =  (47-vdet  C)  exp 


-  \  (i-l)TC_1(i-s) 


,-1 


=  [4"“  ( l-r~-  >>“)]" 

>  exp  |  -  ^ - j  cos?i)2  +  (a2->/2oi  sin^i)2 


cos<  cos$2) 

■  2 x ( a ! - *'FT\  cos<t1)(a4-/^  sin$2) 

+2 x(a2-v^rT1  Sin^Ho^-*^  COS<f>2) 

-2r(a:-  ^/Fc~1  sint  i)(a4-r2cT  sin<t2) 
+(a;.-*/2T:  cos:2)2  +  (a4 -»/2p^  sine2)2 


(2. l-36b) 
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Since  we  are  interested  in  the  phases  at  the  sample  times,  it  is 
appropriate  to  change  from  rectangular  to  polar  coordinates  using  the 
transformation 


Ri  >  0 


|  | 

R2  i  0 

I *2-*2  i 


<  IT 


If  this  transformation  is  employed,  we  obtain  the  joint  pdf  of  envelopes  and 
phases  given  by 


PE P ( R 1  ’  R:*  ^  )  =  RiR2Px(Ricos^i,  R-  sin«i,  R2  cos«2,  R2  sin«2) 


R-  R' 

- — -  exp  [(Rl+R2)/2  +  Pi  +  p 2 

40(1-'-2)  1-p2  L 


-uRiR:  cos( 4 1- ^2 ^ )  -XR 1  cos( ^i-v) 

-Y  R2  COS (  * 2“W)  -  2l‘/C1P2  cos(  O-y+oJ  J 
in  which  we  use  the  notations 

v2  =  O+a2,  i  =  tan~*(yr) 

XJ  =  2C i+2’pzP2  -  4 p 1  p 2  cos( 0i-4>2+^) 

tan  v  =  [>  2c f  sinvi-u^fc^  sin(P2-0]/|}r?c’i  cos<+ i-l»/2F2 


(2.1-37a) 


(2.1-37b) 

(2.1-37c) 
( 2 . 1-37  d ) 
cos(c:-  ) 


(2.1- 37  e ) 
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Y  =  2f;  +  2w‘ p i  -  4u^pip2  cos(^i-4'2+^) 


(2.1- 37  f ) 


tan  w  =  [*Tp2  sin pi  sin(4i +  £)]/[ /2p^  cos$2-v>/2pT  cosUi+S)]. 


( 2 . 1-37  g ) 


As  a  check,  we  note  that  for  y  =  0  (no  noise  correlation),  the  pdf  in 
(2.1-37)  reduces  to  the  product  of  the  pdfs  for  each  sample  time: 


PrP(Ri »  R  >  »  v  i  u=0) 


=  p(Rj » «i  )  p(R.,  «-2  ) 


( 2 . 1- 38a )  i? 


where  each  pdf  is  of  the  form 


P(R,0  =  -jjr  exp  -  -y  -  [  +  R/2 T  cos(4 -4 ) 


(2 . 1  -  38b ) 


as  discussed  in  connection  with  equation  (2.1-30). 

Having  expressed  the  joint  pdf  for  the  envelopes  and  phases  in  (2.1-37), 
the  general  procedure  is  to  integrate  over  the  envelope  variables  to  give  a 
phase-only  pdf, 


P<(-i  .O  =  dRi 


dR  2  P^p(R] »  R; ,  4’i  »  4'2  ) , 


(2.1-39) 


and  then  to  find  the  pdf  for  the  differential  phase  A4  =  4]  -  42.  Various 
approaches  have  been  used  for  carrying  out  the  integrations  shown  in  (2.1-39). 
Due  to  the  complexity  of  the  expression,  the  relative  virtue  of  any  particular 
approach  lies  in  the  feasibility  of  computing  the  phase  pdf’s  obtained  by  the 
approach.  Here  we  will  summarize  some  approaches  that  have  been  tried,  and 
present  the  corresponding  results. 


40 


) 


2. 1.5.1  Fourier  Series 


What  we  call  the  "Fourier  Series  approach"  is  one  used  by  Middleton 
[4,  chapter  9].  The  procedure  is  to  utilize  the  Fourier  series  expansion 


.a  cosb 


=  X  £k  coskb* 

k=0 


in  which 


1,  k  =  0 

2,  k  >  0 


(2.1-40a) 


(2. l-40b) 


and  Ik( . )  is  the  modified  Bessel  function  of  the  first  kind,  of  order  k. 

Using  this  expansion  on  the  cos(4m -v)  =  cos($2  +£<i>-v)  and  cos(42 -w)  terms  in  the 
exponent  of  (2.1-37b),  and  integrating  the  resulting  product  of  series  over  4- 
(a  2r  interval)  yields  a  single  series: 


XR  j  YR: 

d:  exp  -  cos(42  +^4~v)  + - 1  cos(42  -w) 

I--2  1-u2 


■*ri  £^S) 


(2.1-41) 


Next  in  the  procedure  is  to  expand  the  Ri  R2  term  in  the  exponent  to  obtain 


exp - -  cos(£4-5)i 

i-c- 


■  I 


r0 


or  cc 

r— •  r->  u  Rj  R2 

^  ^  1  2(  l-lJ 

£=0  m=0  v 


2m  +  £ 


COSt (A4 -r  )  , 


(2.1-42) 
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in  which  a  series  expansion  for  the  Bessel  function  has  been  used.  At 
this  point,  the  integrations  over  Ri  and  R2  can  be  carried  out.  A  synthesis 
of  formulas  6.643.2  and  9.220  in  [5]  provides  the  integration  formula  needed: 


f  dR  R1+2m+i  e-aR2/2  1  k(s  R^2) 

=  J  dx  xm+i/2  e”ax  I k(2S'/x)  *  2m+*/2 


■rm  +  1  +  (k+t)/21  ,k  -m-l-(k+*)/2  „m+ty2 

k!  -  a  * 


*  iFi[m+l+(k+t)/2;  k+1;  £“/»], 


(2.1-43) 


where  :(•)  is  the  gamma  function  and  iFi(a;b;c)  is  the  confluent  hypergeometric 
function.  After  using  this  formula  twice,  the  pdf  for  44  is  a  triple  infinite 
series  of  the  Fourier  type: 


00  oc  oc 

p..(x)  =  Eli  K(k,*-,m)  cosk(44 -v+w)  cos*.(M-0» 
k=0  £=0  m=0 


where  the  coefficients  are  given  by 


K(k,£,m)  =  -  exp 

2r(l-u2) 


c2+  c2-  v'fp1c2  cos(A<t+5) 

—  — 


k  “i 


m !  ( m+i  )  I  k  1  k  I 


r-  [m+l+(k+*)/2] 


(2. l-44a) 


X 


.  2m+£ 

k* 


>  iFi  [m+l+(k+£ )/2;  k+1;  X2/2(l-'~2)] 

>•  iFi  [m+l+(k+<  )/2;  k+1;  Y2/2(1-l: )]. 

42 


•%.  *v 


(2. l-44b) 
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As  a  special  case,  note  that  for  no  correlation  (y=0),  the  pdf  reduces  to  the 
m=*=0  case  of  a  single  series: 

pjx)  *  ~Z -  T*  £k  (P1P2 )k/2  cos(a$-a») 

k=0  (kl) 

x  ]F1(l+k/2;k+l;Pi)1F1(l+k/2;k+l;p2).  (2.1-45) 

This  special  case  has  been  obtained  by  Mizuno  et  al  [6],  and  is  recognizable  as 
the  self-convolution  of  Middleton's  series  form  [4]  for  the  single  phase  pdf 
(2.1-30),  often  cited  in  the  communications  literature. 

2. 1.5. 2  2tt  Modularity  Issues 

The  procedure  used  in  deriving  the  differential  phase  pdf  (2.1-44)  involved 
a  subtle  but  very  significant  assumption.  In  (2.1-41),  integration  over  $2  was 
assumed  to  be  performed  over  a  2*  interval.  As  presented,  this  integration 
preceded  integration  over  the  envelope  variables,  but  the  order  of  integration 
is  not  an  issue.  What  is  important  is  that  this  step  is  not  strictly  correct 
unless  the  resulting  differential  phase  is  interpreted  as  a  modulo  2ir  quantity. 

To  prove  the  preceding  assertion,  consider  the  following  reasoning.  First, 
we  note  that  the  mathematical  expression  for  the  joint  pdf  P^Ui  ,*2)  is 
periodic,  that  is, 

f^ ($ii  2n jt,  $2  i.  2mrr)  =  f^(<fi»<t2)»  (2.1-46) 

where  f^  is  the  expression.  Next,  we  state  that  the  joint  pdf  consists  of 
the  mathematical  expression  plus  the  restriction  of  4>i  and  $2  to  some  principal 
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V?  ’ 

L*^3  i 

I 


interval  (otherwise,  the  probability  "mass"  under  the  pdf  would  be  infinite). 
Since  for  no  noise,  -»•  C;  and  $2  ^  $2  >  it  is  reasonable  to  state  that 


p.  («*!  ,  *2)  = 


f^i  ,4>2),  1 4>1-  4>  l  |<*  and  I  _  4>2  1  <  * 


0,  otherwise. 


Now,  the  formal  procedure  for  obtaining  the  pdf  of  a$>  =  <J>i  -  <i>?  consists 
of  performing  the  integral 


•U 


P  (y)  =  /  dx  pjx+y,  x)  ,  |a4-ac|  <  2 r, 
l 

where  the  limits  of  integration  are,  by  (2.1-47), 


U  =  min[:2+r,  ;:+~  -  (a$-a<*)] 


and 


L  =  max[::-r,  ::-7  +  (Ac-A4)], 


This  conventional  procedure  yields  a  pdf  pA(J.(')  which  is  nonzero  on  a  47  interval. 
For  example,  for  noise  only  and  no  correlation,  the  "actual"  differential  phase 
pdf  is  triangular: 

(2-  -  |4f  |  )/4 r2  ,  |At|  <  27; 

(2.1-49: 

0,  otherwise. 


p  (y;A=-=0)  = 


But  we  can  also  speak  of  a  "modulo  2*"  differential  phase,  defined  as 

=  (iv  -  A:) mod  2r  +  Ac 


(2 .  l-50a ; 


L?  ,  lAi  -  A  t  I  <  2’T 
(  L\  -27,  27  <  A 4—  Ac  <  4- 

■.1  +27  ,  -47  <  A4  -  A:  <  -2T  . 


(2.1- 50b ; 


(2.1-47)  :-r4  . 


'  V,  * 


( 2 . l-48a ) 


( 2 . l-48b ) 


7 

f.V 


Wf 


(2.1-48c) 


d 


ft  a  ** 

£ 
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The  pdf  for  this  modulo  2t  differential  phase  is([7],  [8])  the  aliased  "actual" 
pdf: 


P  ,(y)  =  1  kr"° 


Y  (y+2k7r) ,  |y-a*|  <  V  ; 


0  ,  otherwise. 


(2.1-500 


Again  using  the  noise-only  example,  modulo  2t  the  pdf  for  the  differential  phase 
is  uniform: 


p,  (y;A=-=o)  =  27  ,  M  <  t  . 


(2.1-51} 


With  these  considerations  in  mind,  we  can  show  the  relationship  of  the 
modulo  2 t  differential  phase  pdf  to  the  integral  form  of  the  "actual"  Li  pdf. 
Using  (2.1-50c)  and  (2.1-48),  we  find  that 


P  (y)  = 


dx  f,(x+y-2r,x) 


ru  o  ru  i 

+  \  dx  ^(x+y,  x)  +  /  dx  f,.  (x+y+2~  ,  x) , 


(2 .  l-52a) 


where 


Jk  =  min[c:+r,  ?z+r  -(y-Ac+2kiT)] 


L,  =  max[;:  -r,  c2-r  +  (-la  -y-2k~) ] . 


(2.1 -52b ) 


(2.1-520 
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For  y(  =  tt)  <  A4  we  have  U_^  <  L_1  ,  so  that  the  first  integral  in  (2.1-52a)  is 


zero,  and 


Uq  =  *2+*  ,  Lj  =  4>2-'it 


Lq  =  Ui  =  <t>2~  ^ 


(2.1-53a) 


(2. l-53b) 


Since  f,  ( • )  is  periodic,  we  therefore  have 


p  (y)  =  /  dx  f4  (x+y,  x) ,  |y-A<?  |  <  -  , 


with  the  same  result  for  y  >  u. 


(2.1-54) 


Our  conclusion  therefore  is  that  integration  over  a  2^  interval  of  $2  in 
the  derivation  of  the  differential  phase  pdf  produces  a  result  that  pertains  not 
to  the  actual  differential  phase  observed  at  the  digital  FM  demodulator  output, 
but  to  a  modulo  2r  version  of  it. 

The  question  is,  which  version  of  the  differential  phase  is  appropriate 
for  analyzing  digital  FM  performance?  Not  using  2t  modularity  greatly  complicates 
the  analysis  and  computation  of  the  pdf  (see  [9]  and  [10],  for  example).  However, 
it  is  incorrect  to  say  that  the  voltage  at  the  output  of  the  receiver  is  pro¬ 
portional  to  a  modulo  differential  phase,  because  there  is  no  mechanism  in 
the  discriminator  or  in  the  integrate-and-dump  filter  which  would  induce  this 
modularity.  For  high  SNR,  the  choice  is  somewhat  arbitrary,  since  the  unaliased 
pdf  is  negligible  for  j a^ -Ac  j  >  tt.  On  this  account,  we  shall  use  modulo  2r 
expressions  unless  stated  otherwise. 


V* 
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2. 1.5. 3  Characteristic  Function  Approach 

Using  a  modified  characteristic  function  approach,  Pawula,  Rice,  and 
Roberts  [8]  have  developed  expressions  for  the  modulo  2*  differential  phase 
distributions.  The  pertinent  results  are  the  following: 


Vy)  *  Hr 


e“^x^  cos  x 


-V2  [1  -  (r  cosy  +  A  sin  y)  cosx]2 


■  E ( x ) +2 


U-W(  rcosA<t  + Asinie 


(2.1-55a) 


where 


r / „ \  _  U-V  sinx  -  W  cos(A$-y)cos  x 
'  l-(r  cosy  +  Xsiny)cosx 


U  =  (m  +  £>-)/2 


V  *  ( 


w  =  vfH:  =  4l2-V2 


(2.1- 55b ) 


(2.1-55c) 
(2.1 - 55d ) 


(2. l-55e) 


P  r  { •»  -, 


F('-:)  -  F(v  i)  +  1,  vi  <  Ac  <  ^2 ; 


F ( C  2  )  -  F(v  1  )  ,  C  l  >  Ai  or  <1’2  <  Ac  ; 


(2. l-56a ) 


where 


H2  -E(x) 


W  sin(Ai  -’^) 


r  si nu  -  x  cos ^ 


4 t,  U-y  sinx  -  W  cos(A(f  -<^)cosx  1  -  (r  cos^  +  Asin^,)  cosx 


A  derivation  of  (2.1-55)  is  given  in  Appendix  A  which  does  not  use  the 
characteristic  function  approach. 


.  ( 2 .  i  - 
- 
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2.1.6  FM  Noise  Clicks 

As  the  signal  plus  noise  waveform  is  processed  by  the  digital  FM  receiver, 
the  total  phase  derivative  is  extracted  by  the  FM  discriminator.  We  have  seen 
that  this  derivative  can  be  expressed  by  the  sum  of  signal  ($)  and 
noise  (n)  terms.  Recall  that  the  phase  noise  is  expressible  as 


n(t)  =  tan" 


v2(t) 

>j(t)  +  /2p(t) 


(2.1-57) 


Now,  when  >c ( t )  (related  to  the  instantaneous  amplitude  of  the  signal) 
is  -small,  the  probability  that  the  denominator  in  (2.1.57)  can  change 
from  a  positive  value  to  a  negative  value  becomes  significant.  The  impact 
of  this  factor  is  illustrated  in  Figure  2.1-6.  In  part  (a)  of  the  figure, 

‘T  is  relatively  large,  so  that  the  example  random  phase  trajectory  is  confined 
to  the  right-hand  side  of  the  origin.  Part  (b)  shows  that,  for  the  same  sequence 
of  values  of  the  noises  v^t)  and  v2(t),  when  is  small,  an  encirclement  of 
the  origin  occurs. 

A  complete  encirclement  of  the  origin  of  course  results  in  a  rapid  2n  increase 
in  r  ,  or  impulsive  value  of  the  derivative  n(t)  -  heard  as  a  "click"  in  FM  re¬ 
ceivers.  Once  the  denominator  of  (2.1-57)  becomes  negative  and  also  the  numerator 
changes  sign,  it  is  highly  likely  that  a  complete  encirclement  will  follow. 
Therefore,  the  expected  number  of  encirclements  can  be  computed  as  the  average 
number  of  zeros  of  v2(t),  given  that  Vj(t)  +  */5p  is  negative. 

Assuming  that  a  click  or  encirclement  always  results  when  ( n |  exceeds  r  , 
we  can  estimate  the  probable  number  of  clicks  per  unit  time  as  follows.  With 
the  help  of  Figure  2.1-7,  we  understand  that  a  positive  click  will  occur  in  the 


numerator 
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interval  (0,  dt)  if  the  phase  noise  angle  n  is  in  the  second  quadrant, 
and  if  it  is  increasing  fast  enough  to  cross  into  the  third  quadrant  during 
the  interval.  In  terms  of  the  quantities  in  (2.1-57)  and  their  derivatives,  the 
probability  of  a  positive  click  in  (0,  dt)  can  be  expressed  as 

P+  =  Pr{v!  <  -  /2c ,  v2  >  0,  v2  +  v2dt  <  0  ) 


f  *1 


-v-dt 


dv2  p(vj  ,V2  ,\>2  )  . 


(2.1-58) 


Similarly,  a  negative  click  will  result  when  the  phase  noise  goes  from 
the  third  to  the  second  quadrant,  with  probability 


P_  =  Pr. <  0,  +  v2dt  >  0> 


d  •  l  J  I  d.’:  p(v1,v2,v2). 

J0  -v- dt 


(2.1-59) 


In  order  to  calculate  these  probabi 1 i ties ,  it  is  necessary  to  develop  the  joint 
pdf  in  (2.1-58)  and  (2.1-59). 

2. 1.6.1  Calculation  of  Click  Rates 

Recall  that  the  total  normalized  waveform  is 


[/2c  (t)  +  v:(t)]  cos[>ot  +  *(t)] 


(t)  sin(>Gt  +  <t(t)] 


=  /2c  (t)  cos[-ot  +  ( t ) ] 


+  •  (t)  cos  a;-  pt  -  v  (t)  sin^0t. 


( 2 . 1 -60a ) 


(2.1-60b) 


with  (.c,-.  )  being  the  conventional  (normalized)  bandpass  noise  components  with 


W  *.  V  o 
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(assuming  a  symmetrical  passband) 


E{vc(t)  vc(  t+T) }  =  E{vs(t)v$(t+T)  =  r(t) 


E{vc(t)  vs(t+T)}  =  -  E  {v$(t)  vc(t+i) }  =  0. 


^(t)  =  vc(t)  cosi(t)  +  v$(t)  s i n <f>( t ) 
-2(t)  =  -vc(t)  sin:(t)  +  v$(t)  cosi(t) 


and  the  noise  correlation  functions  are 


E -^jlthgU+T)}  =  r( t )  sin[c(t)-y(t+i)] . 
E ■-■^1(t)--1(t+T):  =  r( t )  cos[: ( t)-:( t+x ) ] 

=  E  -  >•  2 ( t )  ( t+T ) } . 


From  these  expressions  we  can  deduce  that  the  nonzero  moments  are 


=  Elv2  (t)  ;  =  1  =  Eivfvt)} 


5  E ! v A ( t) !: ( t) }  =  -<(t) 


-  2  E{ v:  ( t) )  =  -  — -  r(x  )  1T=0  +  [*(t )]Z. 


Since  vi,  v:,  and  v:  are  Gaussian  random  variables,  their  joint  pdf  is 


/  *  \  _  1  -vt/2  /  •  \ 

Pi ( • l »- -  * )  j—-  e  P2(vi »v2 )  * 

V  C-7 


(2.1 -61a ) 


(2 . 1 -6 1 b ) 


(2.1-62a) 


( 2 . l-62b) 


(2.1-63a) 


(2.1-63b) 


(2.1-64a) 

'A*. 

(2.1  -  64b )  & 


(2.1-64c) 


( 2 . l-65a ) 
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P p ( ’•'  1  » V 2  )  =  - T==exp - —  V  j 

2  2tto2/i-^  I  2 ( 1-C2 )  [ 


-  2g(— )  +  (-)"  . 


( 2 . 1 -65b ) 


With  this  joint  pdf  we  can  now  calculate  the  click  probabilities  P+  and  P_ 
First  we  note  that 

"Vndt  r  T 


)Ql 

dv:  P1(vi»v2,v2)  =  £  -  Q(~v2dt)  Pp(vi >v2 ) 


/  v2dt  \ 

'  2  er,(^f-  )  P2(v‘^2) 


(2.1-66a) 


(2.1 -66b ) 


/  d'-:Pi  (■> :  »v:  »‘*i )  =  i  -  Q(*:dt)  Pp(v}  ,v; 
-•^dt  1  L  J 


-  2  erf  -V  P2(v^°2)  • 


(2.1-66c) 


(2.1  - 66d ) 


Anticipating  that  we  will  make  the  value  of  dt  very  small,  a  suitable  approxi¬ 
mation  is  the  first  order  MacLaurin  series 


Q(x)  =  \  -  x  •  ^ 


2 

-xV2 


(2.1-67) 


Substitution  of  this  approximation  in  (2.1-66)  and  then  in  (2.1-58)  and 
(2.1-59)  gives 


r° 

\  =  /  d- :  /  dv2  (-v:,dt) 

—  ex  f . 


-(v2dt)V2 


Pp(vl  »  v2 ) 


(2 . l-66a ) 


-(v2dt)  /2 

d-.  (’  ;dt)  '  Pp(v)  >-2 ) 

*  2  r 


(2.  l-68b) 


—  cx 
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} 


are  significant.  In  (2.1-71)  also,  N,  the  net  number  of  opposing  clicks, 
is  assumed  to  have  the  Poisson  distribution: 

Pr(N  =  n)  =  exp{-|Nc|)  •  |NC (n/n!  .  (2.1 

Figure  2.1-8  illustrates  the  effect  of  opposing  clicks  on  the  differential 
phase  probability  distribution. 
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2.2  ERROR  PROBABILITY  FORMULATION 

In  this  subsection,  we  discuss  the  bit  decision  procedures  assumed  to  be 
implemented  by  the  FH/CPFSK  receiver,  and  derive  the  basic  expressions  needed 
to  evaluate  uncoded  bit  error  probability  (BER)  in  consideration  of  thermal  noise, 
clicks,  jamming  noise,  and  intersymbol  interference. 

2.2.1  Decision  Processing 

In  Section  2.1  we  have  discussed  how  the  limiter-discriminator  receiver 
with  integrate-and-dump  filtering  develops  samples  of  the  differential  phase. 

We  denote  those  samples  by{i?  ^  ,i=l ,2,. . . ,L >,  where  L  is  number  of  hops  on  which 
the  data  is  repeated  for  possible  diversity  improvement  against  the  jamming. 

If  L>1,  the  T-second  period  over  which  each  sample  is  developed  is  related  to 
the  bit  period  by  T  =  Tb/L.  In  concept,  we  can  diagram  the  decision  processing 
to  be  evaluated  in  this  report  as  shown  in  Figure  2.2- 1(a):  the  samples  are 
first  processed  by  a  pre-combining  function  g(  •)  to  produce  modified  samples 
zi  ■  ,  where 

z  =g(-i)-  (2.2-1) 

For  example,  g(-)  may  represent  analog-to-digi tal  conversion  with  a  specific 
number  of  quantization  levels;  if  the  number  of  levels  is  just  two,  then  the 
function  g(. )  implements  the  "hard  decision"  given  by 


d  *  0 


(a)  Conceptual  Form  of  Decision  Processing 


j-th  sample  on  hop  i 


(b)  Details  for  Slow-hopping  System 


FIGURE  2.2-1 


RECEIVER  DECISION  PROCESSING 
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The  modified  samples  are  then  summed  to  produce  a  single  decision  statistic 


z,  that  is, 


2  -  Z  2t =  Z  s(2V’ 


so  that  the  final  bit  decision  d  is  made  by  examining  the  sign  of  z: 


1,  z  5  0 


0,  z  *  0. 


In  practice,  since  the  data  rate  is  assumed  to  be  higher  than  the 
hopping  rate,  the  L  "chips"  constituting  the  repetitions  of  each  "data 
bit"  are  not  transmitted  one  after  another,  as  implied  by  Figure  2.2-l(a). 
Instead,  they  are  transmitted  as  members  of  sequences  of  Q  chips  per 
hop,  and  the  sequences  are  repeated  over  L  different  hops  with  the  idea  that 
some  of  the  hops  will  evade  the  jamming.  Thus,  some  prearranged  formatting 
is  necessary  so  that  the  receiver  can  retrieve  the  L  "pieces"  of  a  given  bit. 
Figure  2 . 2- 1  ( b )  illustrates  the  processing  necessary  to  collect  a  particular 
data  bit's  L  "pieces"  from  among  the  set  of  LQ  modified  samples  which  are 
received  over  the  L  hops.  For  our  analyses,  it  is  not  necessary  to  describe 
the  receiver  decision  processing  any  further,  except  to  note  that  scrambling 
or  re-ordering  of  the  Q-bit  sequence  from  hop  to  hop  is  an  option;  if  such  a 
re-ordering  is  done,  then  the  intersymbol  interference  differs  from  hop  to  hop, 
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and  this  affects  the  analysis,  as  we  shall  mention  below  at  the  appropriate 
place. 

2.2.2  Conditional  Error  Probabilities 

Each  of  a  data  bit's  L  differential  phase  samples  M.  is  a  random  variable 
whose  distribution  is  parametric  in  the  values  of  the  following  quantities 
during  the  sample  interval:  the  CNR  p^,  representing  the  thermal  noise  and 
the  presence  or  absence  of  jamming;  and  the  nominal  differential  phase  a<*.  and 
the  CNR-related  parameters  U-j ,  V-j ,  and  ,  representing  the  data  pattern- 
dependent  intersymbol  interference  effects.  Thus  the  distribution  of  the 
modified  sample  zi  is  conditioned  on  these  quantities  as  well: 


P2.(;H-(i),  sp  =  pz.(^k(i),  ui,vi,w.,A^i), 


(2.2-5) 


where  the  vector  is  a  shorthand  notation  for  the  set  of  data-dependent 
parameters  listed. 

The  symmetries  of  the  data-dependent  parameters  are  such  that 


xly}  =  Pr{A$.j  >  0[p^>  x  Oy  }, 


where  "xly"  stands  for  a  data  pattern  with  a  data  value  of  1  at  the  sample 


(2.2-6; 


time,  and  "xOy"is  the  complementary  pattern.  If  the  pre-combining  processing 
function  g(-)  is  perfectly  symmetric  (odd),  then  it  is  also  true  that 


Pr{z-  <  0 i xly}  =  Pr{zi  >  0|x0y 


(2.2-7) 


and  for  L=1  we  may  restrict  our  attention  to  the  xly  patterns  in  calculating 
the  BER. 


/■A-'vVvVvvV'i-v'-V 

A  -  ■  Mu', 


For  L>1,  the  randomness  of  the  { z ^ )  is  due  to  noise  power  at  different 
portions  of  the  RF  spectrum,  so  they  are  statistically  independent,  with  their 
joint  pdf  given  by 


IT  Pz  .(“i  *  i-j)  • 

i  =  l  1 


(2.2-8) 


The  parameter  sets  {c^,  6.}  are  in  general  different  for  each  hop,  with  the 
L  values  of  independent  of  each  other  since  it  is  assumed  that  the  event 
of  jamming  on  one  hop  is  independent  of  the  event  of  jamming  on  another  hop. 

The  data  dependent  parameter  values,  however,  are  related,  since  the  data 
sample  values  are  the  same.  There  are  two  cases  we  shall  consider:  (1)  the 
transmitted  data  sequences  on  the  L  hops  are  identical,  and  (2)  the  transmitted 
data  sequences  are  scrambled.  For  case  (1),  the  data-dependent  parameters  are 
identical  i  ,  all  i).  For  case  (2),  the  individual  are  independently 

selected  from  those  corresponding  to  the  four  possible  data  patterns  represented 
by  xly  ,  or  xOy  . 

In  consideration  of  the  symmetries  we  have  discussed  with  respect  to  the 
distribution  of  individual  hop  samples,  it  is  possible  to  state  the  symmetry 
which  applies  to  the  conditional  probability  of  error  as  follows.  The 
distribution  of  the  sum  decision  statistic  is  such  that  when  x  =  0 


Pr{z  >  Oj;^,  j32;...;p{L)  _6 L } , 

=  Pr(z  <  0|p(1\  ?i;p(2),  £,_> ,  (2.2-9) 

where  denotes  the  data-dependent  parameter  values  arising  from  the  comple¬ 
ment  of  the  pattern  which  produces  .  Therefore,  the  conditional  probabilities 
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of  error  have  the  symmetric  relationship 
P(e|p(1),  x1ly1;...;P^L),xLlyL) 


a 

& 

8 


P(e]p^  .XjOyjj. . .  ;p^  >xL0yL) , 


(2.2-10) 


when  x  =  0,  which  holds  for  a  noise  spectrum  symmetric  about  the  carrier  frequency. 

2.2.3  Unconditional  Error  Probability 

The  unconditional  error  probability  is  found  by  averaging  the  conditional 
error  probability  with  respect  to  the  jamming  events  and  the  intersymbol  inter¬ 
ference  patterns.  The  required  pdfs  for  averaging  are  assumed  to  be  discrete¬ 
valued. 

The  averaging  over  the  jamming  events  is  most  easily  done  at  the  per-hop 
level.  That  is,  in  the  calculation  of  the  error,  we  use  the  hop  pdf's 


pz  (all,)  =  Y.  Pr{°^  =  R  }  pz  ij) 


(2.2-lla) 


(1-Y)  Pz  (a|pN,  6.)  +  YPZ  (a|pT,  ) , 


(2.2-llb) 


where  the  unjammed  SNR  and  the  jammed  SNR  py  are  related  to  the  signal 
bit  energy  and  noise  spectral  densities  by 


(2.2-12a) 


I  .  __b _ 

L  Nn+Nn , 


1  .  b 

L  Nn+N Jy 


(2.2- 12b) 
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In  ( 2 . 2- 12b ) ,  it  is  assumed  that  the  actual  jamming  spectral  density  Nqj  is 
related  to  the  average  jamming  spectral  density  Nj  by  Nqj  =  N j/y -  This  is 
turn  assumes  that  the  fraction  y  of  the  hopping  band  is  jammed  with  noise 
density  Nqj;  the  fraction  y  also  is  the  probability  of  jamming,  as  used  in 
(2.2-llb) . 

Assuming  that  the  partially-averaged  pdf's  given  in  (2.2-11)  are  used, 
the  conditional  probabilities  of  concern  are  of  the  form 

P(e | Xj lyj , . . . ,x^ly^) . 

The  complementary  patterns  need  not  be  considered  in  view  of  the  symmetries 
discussed  above,  assuming  that  1 * s  and  0's  are  equiprobable.  Therefore, 
the  general  BER  is  the  average 


(2.2-13) 


p(e)  =  Y.  (i  )  P(e‘x11yi>---»xLiyL)> 


(2.2-14) 


if  the  data  sequences  are  scrambled  from  hop  to  hop.  If  identical  data 
sequences  are  used,  then  =  x  and  y^  =  y,  and 

P(e)  =  P( e i 010)  +  P ( e 1  Oil)  +  P(e| 110)  +  P(e|lll)}. 

As  a  practical  matter,  if  the  sequences  are  scrambled,  it  is  more 
convenient  to  use  the  marginal  (averaged)  sample  pdf’s 

Pz(y)  =  Pz .[c ! 6(010)]  +  PZT»|6(011)] 

+  p_  [,|i(110)]  +  p  [a  |  6_(  1 1 1 )  ]  >  , 

L  i  4  i 

and  to  straightforwardly  calculate 

P(e)  =  Pr{^  <  0|1>. 


(2.2-15) 


(2.2-16) 


(2.2-17)  ; 
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2.2.4  Treatment  of  Cl icks 

The  single-hop  pdf  given  by  (2.2-11)  may  be  expanded  in  terms  of  the 
number  of  FM  noise  clicks,  as  follows: 

P2  (aU)  =  (l-Y)pz.MDN,£)  +Yp2  (a|pT,B) 


(l-'O  Y,  Pr{Nc=nl°N*i)}  f(a+2^n;pN,B) 


+  i  Y  PrtNc=n  |  f  (a+2Trn;Py  ,_e) , 


(2.2-18) 


with  Nc  the  number  of  clicks.  For  convenience  we  consider  only  patterns  with 
>  0,  so  that  the  probable  clicks  of  significance  are  negative,  those  which 
have  the  effect  of  shifting  the  differential  phase  pdf  to  the  left  by  some 
multiple  of  2r .  In  (2.2-18)  we  indicate  an  infinite  sum  with  respect  to  the 
number  of  clicks  in  the  interval  T;  this  can  be  interpreted  as  averaging  the 
pdf  f(,+2rN;.  ,^)  with  respect  to  N.  We  also,  in  using  £_  (without  a  subscript), 
suppose  that  no  scrambling  of  the  data  is  done  from  hop  to  hop;  then  (2.2-18)  is 
the  pdf  for  each  of  the  L  sampled  differential  phases,  conditioned  on  the  data 
pattern. 

2.2.4. 1  Number  of  Significant  Clicks 


It  is  possible  to  write  (2.2-18)  as 


P2  LiL>  =  Y  fn^ 


(2.2-19) 


emphasizing  the  fact  that  this  pdf  is  the  superposition  of  Zv  translations  of 
pdfs  for  given  numbers  of  clicks.  We  then  can  appreciate  that  the  sum  z  also 
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has  a  pdf  which  is  of  this  form,  that  is,  is  a  superposition  of  (overlapping) 
pdf's  for  given  numbers  of  clicks.  When  we  consider  that  each  of  these 
constituent  pdf's  is  nonzero  for  only  a  finite  portion  of  the  2-axis,  then  it 
follows  that  a  certain,  finite  number  are  nonzero  for  positive  z.  This  fact 
allows  us  to  restrict  our  attention  to  a  finite  number  of  clicks  when  computing 
the  probability  of  error. 

The  reasoning  is  as  follows:  the  sum's  pdf  may  be  expressed 
by 


pz ( -  !.£.)  =  Y  fn^  (a+2lTn;D 

n=0 

co 

=  )  g  [a-Uc-+2*n;B],  (2.2-20) 

n=0  n 

in  which  the  gn(  :*)  are  pdf's  centered  at  a  =  0,  so  that  the  constituent  pdf's 
f^(')  are  seen  to  be  centered  at  LA;  -2a n. 

If  the  single-hop  pdf's  are  nonzero  for  a  2r  interval,  then  it  follows 
that  the  gn  are  nonzero  for  2tI  intervals  as  a  result  of  convolution  of  L 
pdf's.  Now,  the  error  probability  can  be  computed  from 

P(e)  =  Pr{z<0) 

=  1  -  Pr(z>0).  (2.2-21) 

Therefore,  in  computing  P(e)  we  need  only  those  constituent  pdf's  which  are 
nonzero  fur  l  -  0,  corresponding  to  click  numbers  n  =  0  to  n  =  nmax,  where 


i 
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nmax  is  deduced  from 


LAs}>  -  2ir  nmax  +  u  <  0 


nmax  >  (a<j>+  ^) . 


(2.2-22a) 


For  example  if  =  tt/2,  then  nmax  is  the  smallest  integer  larger  than  3L/4. 


It  is  further  deduced  that 


[&»♦¥»] 


L-l ,  A* /it  >  ^  , 

L-2,  A$/n  <  , 


( 2 . 2- 22b ) 


for  0  <  A<f>  <  ir. 


Figure  2.2-2  illustrates  for  A$  =  .6u  the  facts  that  for  L=l,  the  P(e) 
can  be  computed  using  just  the  pdf  for  no  clicks;  for  L=2,  using  the  sum  pdf's 
for  zero  and  one  click;  and  for  L=3,  using  zero,  one,  and  two  clicks. 

Knowing  this  basic  rule  on  the  number  of  clicks  that  are  significant 
for  computing  the  error  saves  much  computation,  and  also  permits  simplification 
of  the  analysis  somewhat. 

2. 2. 4. 2  Application  to  Sum  Characteristic  Function 

With  the  preceding  in  mind  we  can  truncate  the  click-indexed  series 
(2.2-19)  and  write  the  single-hop  pdf  as 


™  (1) 

Pz-(a|£)  =  2_  fn  (a+21Tn;6), 


(2.2-23) 
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with  characteristic  function 

(1)  n ^ 


nmax 

(v).  £ 

n=0 


where 


(2.2-24a) 


(2.2-24b) 


Then,  since  the  samples  are  identically  distributed  (given  £) ,  the  characteristic 
function  for  the  sum  of  L  samples  is 


JLk.)  =  [^(1)(v)]L 

r nmax  „  , 

■[£  *„<->*'J2’nvL 

n=0 

=  7  (,  ,  L 

\  £-o»Li  *•••>£ 


i:  =  VL 


L  \  nmax  tk  -jZr-vkJ.^ 

£0  ,Li  , . . .  *£  nmax  /  TT  k 


(2.2-25) 


The  first  several  terms  of  (2.2-25)  are 


3.0  BER  CALCULATIONS  FOR  DIVERSITY  SUM 


>; 

4X5 


In  this  section,  we  develop  theoretical  expressions  and  numerical  results 
for  the  BER  produced  by  the  limiter-discriminator  FH/CPFSK  receiver  using 
straightforward  linear  combini  lg  of  the  L  diversity  samples  of  the  differential 
phase.  Worst-case  partial-band  noise  jamming  (WCPBNJ)  is  assumed,  and  our 
interest  is  to  determine  the  improvement  in  the  BER,  if  any,  that  is  achieved 
by  processing  the  diversity  samples  in  this  manner. 

3.1  L  =  1  BER  IN  WORST-CASE  PARTIAL-BAND  NOISE  JAMMING 

The  baseline  case  for  measuring  improvement  in  performance  is  the  case 
of  no  diversity,  that  is,  L  =  1.  The  BER  for  this  case  is 


Pr ;z  =  Li  <  0  |  xly } 


x.y 


(3.1-1) 


where  as  before  "xly"  denotes  the  different  data  patterns  which  determine  the 
intersymbol  interference  parameters,  and  the  effects  of  jamming  and  noise 
clicks  are  already  included.  The  pdf,  conditioned  on  the  pattern,  is  given  by 


P2(-*  £)  -  (!■')  Pgjn  fU+2-n;oN>6)  +v  p^>n  f(a+2i’n;;:1-,_6)  ( 3 . 1  -2a ) 

n  n 

=  X  fn1J  (^+2rn;s)  s  P^M^.),  ( 3 . 1  - 2b ) 

n 

in  the  notation  of  Section  2.2.4,  with  the  click  probabilities  p^  p  (unjammed) 
and  p.  (jammed)  given  by 

i » n 
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Pr  {  Nc  =  n  j  ,j) }. 


The  conditional  BER  is  conveniently  expressed  as 


P(ejxly)  =  1  -  Pr{z  >  0|xly) 

*  Ad+ir 

=  1-1  d  01  P2(aiB) 

Jo 

=  1  - J  di  f^1* (“;£). 

Note  that  only  the  term  of  the  pdf  corresponding  to  zero  clicks  is  needed 
Expanding  further  yields 


P(e  xly)  =  1  -  (l->)P0jo 


■  L:+t 

d'>  f(a;cN»£) 


'L:  +  Ti 


-  ■»  p 


’  P  1,0  d  l  f(a;cT^} 


=  l  -  (l-o  P 


0,0 


-  •>  p 


1,0 


-  (l->)p0>0  [F( ^c+tt;cn ,£)  -  F(0;pn,6)] 


p.  n  [F(o-  +tt;dt,£:)  -  F(0;ot,£)], 
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where  the  function  F(-)  is  given  in  { 2 .  l-5b) ,  with  >.  =  0. 

3.1.1  Case  of  No  Jamming 

With  no  jamming,  the  BER  is  (3.1-5)  with  the  substitution  of  r  =  0, 
resulting  in 

P(eixly)  =  1  -  P0,0  -  P0,0  F(A*+1T;cN’£)  +  p0,0  F(°iPN^- 

Averaging  over  the  data  patterns  (with  parameters  denoted  by  _6)  gives 

P(e)  =  |  ' P(e ; 010)  +  P(e ! Oil)  +  P(e 1 110)  +  P(e!lll)}. 

3 . 1 . 1 . 1  Pattern-dependent  Quanti ties 

To  illustrate  the  use  of  the  various  parameters  developed  in  Section 
2.1,  we  list  in  full  the  pattern-dependent  quantities: 


(3.1-6) 


(3.1-7) 


Pn  n  =  e 


(3.1-8a) 


where 


zi  I  dt  u(t)v(t)  -  u ( t ) v i 
27  _T  u2(t)  +  v2(t) 


exp  {  -c[u2(t)  +  v2(t)]} 


and  u(t)  =  c1  cos(rt/T) 

v(t)  =  c2  "  c3  cos(2~t/T) 


(3.1 -8b ) 


(3.i-8c) 


(3.1-8d) 


For  computation  it  is  convenient  to  define  the  following  functions: 
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u(t)  =  u(Tt) ,  v ( t )  =  v(Tt) 


(3.1-9a) 


w(t)  =  Tu(Tt)  ,  z(t)  =  Tv(Tt). 


(3. l-9b) 


With  these  functions,  ( 3 . 1 -8b )  becomes 


N  dx  w(x).v(_x  ^uixMxj 

C  2r  4  ( u'(  X ) ) 2  +  (v(x))2 


exp{- o[(u(x))2  +  (v(x))2] }  . 


(3.1-10) 


The  functions  u',  v",  w,  and  z  for  the  several  patterns  are  given  in  Table 


3.1-1. 


Now,  the  parameter  set  includes  the  nominal  differential  phase  , 
and  the  CNR-related  parameters  U,  V,  and  W  defined  in  Table  2.1-3.  For  the 
010  pattern,  these  are 


( 0 i 0 )  =  2  tan'1  |c  _1-|  =  1.2239  radians 
U(010)  =  [c2  +  (c2-c3)2]  =  0.8696s 


V(010)  =  0 


(3.1-11: 


W(010)  =  U(010). 


The  numerical  values  assume  that  h  -  0.7,  WjpT  =  1.0,  and  a  Gaussian  shaped 


I .F.  filter. 


PATTERN-DEPENDENT  FUNCTIONS  FOR  COMPUTING  THE  ERROR  PROBABILITY 
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To  calculate  P(e|010),  Pg  g  is  computed  by  numerical  integration  using 
(3.1-8a)  and  (3.1-10);  the  6_  parameters  are  calculated  using  (3.1-11)  and 
Table  2.1-1,  and  substituted  into  numerical  integrations  of  (2.1-56). 

Oil  Pattern 

For  this  pattern  (3.1-8a)  and  ( 3 . l-8b )  apply,  with  u(t)  and  v(t)  given 


u ( t )  =  c4  sin(rt/2T)  -  c^  sin(3rt/T) 


(3.1 - 12a ) 


v(t)  =  C  6  +  C7  COS  (t  t/T )  -  Cg  cos(2rt/T). 


(3.1-12b) 


The  transformed  versions  of  these  functions  and  their  derivatives,  shown  in 
Table  3.1-1,  are  used  in  (3.1-10)  to  compute  the  click  number  average  Nc. 
The  parameters  needed  for  computing  F(>)  are,  from  Table  2.1-3, 


.:(011)  =  tan 


•1  T  c4+c5 
c6~c7"c8 


=  1.7108  radians 


U(011)  =  .  [(c4+c5)2/2  +  c 2  +  (c6-c8)2]  =  0.7784; 


V(011)  =  .  [2c7(c6-c8)  -  (c4+c5)V2]  =  -  0.0997. 
k(011)  =  V  U2(011)  -  V2(011)  =  0.7719; . 


(3.1-13) 


WV 
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110  Pattern 

The  functions  u(t)  and  v(t)  for  this  pattern  are 


u(t)  =  cos(rt/2T)  +  cos(3rt/2T) 


v(t)  =  c6  -  c-j  cos(r t/T)  -  cg  cos(2n/T), 


with  the  transformed  versions  shown  in  Table  3.1-1.  We  note  that  since, 


( 3 . l-14a ) 


(3.1-1 4b ) 


for  example, 


u( t; 110)  =  u(t+T;011)  =  u(-t-T;011) 


v( t; 110)  =  v(t+T;011)  =  v(-t-T;011), 


( 3 . 1 -15a ) 


(3.1-15b) 


there  is  a  great  deal  of  symmetry  with  the  pattern  Oil.  In  fact,  using 
(3.1-15)  we  see  that 


^(110)  =  —  J  dt  :(-t-T;011)e"Ca  {‘t_T;011) 


=  ~  f  dt  : ( t-T ;011)e-Ca^^t-T’011) 


=  ~  J  dt  ’  (t;011)e'Ta  (t;011)  =  Nc(011).  (3.1-16) 

We  also  note  from  Table  2.1-3  that  L:  ,  U,  and  W  are  the  same  for  pattern 
110  as  for  pattern  Oil.  The  sign  of  V  changes,  but  this  is  of  no  significance 
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when  one  examines  equation  (2.1-56)  under  a  change  of  the  variable  of 
integration  from  x  to  -x.  Therefore,  we  can  conclude  that 


P(e - 110)  =  P(e;011), 


(3.1-17) 


and  avoid  the  necessity  for  computing  P( e | 110)  separately. 


Ill  Pattern 

This  "all  l's"  pattern  has  the  simple  result  that 


T  -  h  /:a° 
*c  -  -  2  e 


since  u(t)  =  aQ  sin(rht/T)  and  v(t)  =  a^  cos(^ht/T).  The  pattern 
parameters  for  calculating  F(-)  are 


(3.1-18) 


U(  111)  =  W(  111)  =  a^:  =  0.6806.- 


V(  111 )  =  0 


(3.1-19) 


.:  ( 1 1 1 )  =  rh  =  2.1991  radians. 


In  view  of  the  symmetries  we  have  noted,  the  P(e)  expression  (3.1-7)  can 


be  modified  to 


P(e)  =  i  •;  P(e  010)  +  2P(e  1 0 1 1 )  +  P ( e  1 1 1 1 ) > . 


(3.1-20) 
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3. 1.1.2  Numerical  Results  for  No  Jamming 

The  average  BER  (3.1-20)  is  computed  using  the  program  given  in  Appendix 
B,  with  the  results  as  shown  in  Figures  3.1-1  to  3.1-3. 

In  Figure  3.1-1,  the  time-bandwidth  product  D  =  WjpT  is  fixed  at  the 
value  of  D  =  1,  while  the  digital  FM  modulation  index,  h,  is  varied.  The 
well-known  property  that  h  =  0.7  is  the  best  value  for  the  additive  Gaussian 
channel  with  high  SNR  is  illustrated  by  the  figure.  Crossovers  in  the  figure 
reveal  that  this  accepted  value  of  h  is  best  for  SNR's  greater  than  4-5  dB; 
for  lower  SNR,  h  =  0.8  is  the  best  value  among  the  values  of  this  parameter 
which  are  plotted. 

In  Figure  3.1-2,  the  modulation  index  is  held  fixed  at  the  value  h  =  0.7, 
while  D  is  varied.  Plotted  against  SNR,  as  in  this  figure,  the  BER  decreases 
uniformly  as  the  filter  bandwidth  increases  (D  increases).  However,  this  is  a 
misleading  portrayal  of  the  dependence  of  the  BER  on  D,  since  the  SNR  is 
given  by 


ST 


nowift 


(3.1-21) 


in  order  for  SNR  to  remain  constant  while  D  varies,  the  bit-energy-to-noise 
density  ratio  must  vary.  A  more  fair,  equal-bit-energy  comparison  is  presented 
in  Figure  3.1-3,  in  which  the  BER  is  plotted  against  E^/Nq.  This  second 
comparison  is  more  in  accord  with  intuition,  for  we  observe  a  tradeoff  between 


77 


0.  00 

FIGURE  3.1 


3.00  6.00  9.00  12.00  15.00 

SIGNAL  TO  NOISE  RATIO  (dB) 

-1  UNJAMMED  CPFS K  BER  VS.  SNR  FOR  D  »  1,  WITH  MODULATION  INDEX 
VARIED 


78 


PROBABILITY  OF  BIT  ERROR 


J.  S.  LEE  ASSOCIATES,  INC. 


the  beneficial  effects  of  increasing  0  (more  signal  energy  recovered)  and 
the  degrading  effects  (more  noise  power  admitted),  with  D  =  1.0  being  optimum 
for  high  E^/Nq,  D  =  0.9  best  for  8  <  E^/Nq  <  12  dB,  and  D  =  0.8  best  for 
Eb/NQ  <  8  dB,  among  the  values  of  D  used  in  the  figure. 

3.1.2  Case  of  Partial-Band  Jamming 


The  computations  for  no  jamming  are  parametric  in  the  SNR,  yielding 


P(e)  =  P(e;cN). 


(3.1-22) 


Inclusion  of  partial-band  jamming  is  performed  by  calculating  the  expression 


P(e)  =  (l->)  P(e;.-N)  +>P(ej=T), 


(3.1-23) 


where 


N+  "j 


(3.1-24) 


is  the  effective  SNR  when  the  signal  is  jammed,  and  y  is  the  fraction  of  the 
hopping  band  which  is  jammed. 

Figure  3.1-4  shows  plots  of  the  jammed  BER  (3.1-23)  as  a  function  of 

Eb/Nj  =  .j,  and  parametric  in  v  ,  for  E^Nq  fixed  at  11.75  dB.  This  value  of 

-5 

Efa/No  gives  a  10  BER  without  jamming,  as  can  be  observed  in  the  figure  as 
Eb/Nj  becomes  large.  For  each  value  of  vf  the  BER  curve  is  "S-shaped,"  with 
a  left  (low  Eb/N j  or  strong  jamming)  asymptote  close  to  the  value  of  P(e  )  =  V2, 
and  with  a  right  asymptote  (high  Eb/Nj  or  weak  jamming)  of  P(e)  =  10 


BIT  ENERGY  TO  JAMMING  DENSITY  RATIO  (dB) 

FIGURE  3.1-4  FH/CPFSK  BER  VS.  E / N j  IN  PARTIAL-BAND  JAMMING,  FOR  ONE  HOP/BIT 
AND  E5/Nq  =  11.75  dB  (GIVING  A  10'5  BER  WITHOUT  JAMMING), 
PARAMETRIC  IN  THE  FRACTION  OF  BAND  JAMMED. 
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The  centers  of  the  "S-curves"  are  different  for  each  value  of  y,  so  that 

when  they  are  superimposed  as  in  the  figure,  it  is  evident  that  a  particular 

value  of  y  produces  the  worst-case  or  highest  BER  value  only  for  a  small 

range  of  E^/Nj.  For  example,  the  value  y  =  0.01  is  the  worst-case,  among  the 

values  displayed,  for  21  dB  <  E ^/N j  <  24  dB,  approximately. 

Note  how  much  greater  the  error  is  with  worst-case  jamming  than  with 

full -band  jamming  (y  =  1).  For  example,  10"4  is  achieved  when  y  =  1  for 

-4 

E^/Nj  =  17  dB;  when  y  is  the  worst-case  value,  a  10  BER  requires  Eb/Nj  = 

34  dB,  or  17  dB  more  signal  power. 

If  we  were  to  superimpose  additional  curves  for  more  closely-spaced 
values,  it  is  easy  to  predict  that  the  appearance  of  Figure  3.1-4  would  become 
essentially  that  of  a  black  solid,  with  its  upper  "edge"  a  plot  of  the  worst-case 


BER  for  E^/Nj  continuously  varied,  or 
PWC  =  max  '’('W- 


(3.1-25) 


Figure  3.1-4  suggests  that  the  slope  of  that  worst-case  BER  curve  would  be  -1, 


that  is, 


p«c  *  or  • 2  dB  *  VNJ  * 35  dB 

b  J 


(3.1-26) 


for  much  of  the  range  of  E^/Nj,  with  the  slope  eventually  increasing  to  zero  as 
Ep/Nj  increases  beyond  35  dB  or  so,  since  the  BER  cannot  be  less  than  the 

_5 

un jammed  value  of  10  . 


1 


If  E^/Nq  is  increased  to  15  dB,  giving  an  unjammed  BER  <<  10  (see 
Figure  3.1-1),  essentially  the  jammed  BER  is  given  by 

P(e)  =  Y  P(e;cT) ,  (3.1-27) 

and  the  superimposed  fixed-Y  curves  are  as  shown  in  Figure  3.1-5.  They  are 
still  "S-curve:;"  as  in  Figure  3.1-4,  but  the  high  Eb/Nj  asymptote  is  much  lower, 
and  the  low  Eb/Nj  asymptotes  are  somewhat  lower.  Now  we  see  that  the  worst- 
case  BER  is  characterized  by  a  -1  slope  throughout  the  range  of  Eb/Nj  shown,  so 
that 

pwc  =  yf]  •  W2“B-  <3-1-28a 

with  the  worst-case  value  of  ,  being  very  closely  predicted  by 


WC 


1.58 


Eb/NJ 


Eb/Nd 


1, 


Eb/NJ 


1.58  =  2  dB 

1.58  =  2  dB. 


(3. 1 - 28b ; 


3.1.3  Simplified  Calculations  for  L  =  1 

It  shall  be  interesting  to  compare  the  exact  BER  results  for  L  =  1, 
presented  above,  with  calculations  based  on  simplifying  assumptions.  The 
expressions  (2.1-55)  and  (2.1-56)  for  the  differential  phase  distribution's 
pdf  and  probability  function  simplify  considerably  for  the  following 
assumptions: 


*V'*y*v 
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FIGURE  3.1-5  FH/CPFSK  BER  VS.  Efa/N j  IN  PARTIAL-BAND  JAMMING  FOR  ONE  HOP/BIT 
AND  Eb/NQ  =  15  dB,  PARAMETRIC  IN  \  THE  FRACTION  OF  BAND  JAMMED 
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A;  =  ~/2  (representative  or  typical  value) 


r  =  >.  =  0  (neglecting  correlations) 


V  =  0 


U  =  W  =  0  a* 


(3.1-29) 


(neglecting  intersymbol  interference 
asymmetries) . 


In  effect  we  are  replacing  the  distribution,  averaged  over  the  pattern- 

dependent  quantities,  with  a  distribution  with  "typical"  values.  In  that 

the  V  =  0  assumption  relates  to  the  "all  one's"  pattern,  it  is  reasonable 

to  interpret  i:  =  r/2  as  arising  from  a  modulation  index  value  of  h  =  0.5, 

2 

so  that  a  =  0.81  ;  exp  ;-r/l6:  is  appropriate.  However,  the  analysis 
approach  used  here  is  not  so  formal  as  to  prevent  an  arbitrary  choice  of  some 
value  for  a  ,  if  it  turns  out  to  model  the  exact  results  well. 

The  jammed  BER  for  the  simplifying  assumptions  becomes 


P(e)  =  (1-  •»)  I  cQe  + 


where 


1  "i 

ci  =  j  e  =  P(e! no  clicks,  i ) , 


aTj  tz 


•.“l*  *.*, 

-A 

7a 


>'  sS- 

v  H  A* 

*. 

,VV  s 


-CqI  r  ~ci  “^il 

!  -  e  UJ  +  y  Cje  1  +  1  -  e  ij,  (3.1-30a)  '■ 


(3.1- 30b )  TT.  V 


,  i  =  0  (unjammed) 


T  ,  i  =  1  (jammed) 


-.v  v 


(3.1-30c) 
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By  coincidence,  the  probability  of  zero  clicks  turns  out  to  be  e  '  as 
shown;  that  is,  if  the  phase  modulation  is  assumed  to  be  such  that  *  =  ^h/T 
v/  21. 

O 

For  no  jamming  (r  =  0),  the  simplified  BER  expression  using  a  =  0.81 
gives  the  following  values: 


N 


10  dB 

11.25733  dB 
15  dB 


BER 

1.5176  (-4) 
1.0000  (-5) 
1.8781  (-12)  . 


Calculation  of  the  worst-case  partial-band  jamming  BER  produces  the  following 
asymptotic  results  (higher  Efc/Nj): 

P. 


N 


11.25733  dB 


_WC 
0.23/(Eb/Nj) 


WC 

1.7/(Eb/Nj) 


15  dB 


0.22/(Eb/N0)  1.4/(Eb/Nj). 


(3.1-31a) 

(3.1-31b) 


Note  that  this  approximation  gives  a  higher  jammed  error  than  the  exact 
expression  [see  (3.1-28)],  and  a  lower  unjammed  error.  Thus,  "tinkering" 

P 

with  the  vlaue  of  a  (to  scale  the  SNR)  will  not  accomplish  agreement  between 
approximate  and  exact  answers  for  both  unjammed  and  jammed  conditions. 
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Another  approximation  approach  is  to  fit  a  simplified  curve  to  the 
exact  unjammed  P(e)  shown  in,  for  example,  Figure  3.1-1.  For  h  =  0.7  and 
D  =  1.0,  we  have  the  approximation 


-0.717c 


P(e)  =  0.394  e 


which  is  fitted  precisely  at  o  =  0  dB  and  10  dB;  it  gives  a  10"°  BER  for 


(3.1-32) 


11.7  dB.  When  this  form  is  used  in  the  jammed  BER  equation,  for  high  Eb/Ng 
(negligible  thermal  noise),  the  worst-case  BER  obtained  is 


Vnj 


(3. l-33a ) 


WC  -  Eb/Nj 


(3.1-33b) 


These  results  compare  well  with  the  exact  case  given  in  (3.1-28). 
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3.2  8ER  COMPUTATIONS  FOR  L  ■  2 

The  calculations  of  the  BER  for  L  =  1  showed  that  worst-case  partial-band 
jamming  severely  degrades  the  FH/CPFSK  performance.  We  now  begin  to  determine 
whether  linear  diversity  combining  will  produce  an  improvement. 

3.2.1  Derivation  of  BER  Expressions 

The  two  differential  phase  samples  from  two  hops  for  a  particular  bit  are 
statistically  independent,  due  to  the  separation  in  time  and  frequency  of  the 
noises  added  to  the  signal  on  the  two  hops.  Therefore,  the  joint  pdf  for  the 
two  samples  is,  for  a  particular  data  pattern, 

Pz,7  (x.yii)  =  PU)(x|£)  p(1)(y|£),  (3.2-1 

Z1 z2 

using  the  notation  of  (3.1-2). 

Restricting  our  attention  to  patterns  with  a  bit  value  of  1,  the 
conditional  probability  of  error  is  given  by 

P(e  £)  =  Pnzj+zg  <  Oje  }  ,  (3.2-2 

and  the  probability  of  error  is  represented  by  the  area  noted  in  Figure  3.2-1. 
Since  the  joint  pdf  (3.2-1)  is  symmetric  about  the  line  z 1  =  z^,  the  conditional 
BER,  with  the  help  of  the  figure,  is  seen  to  be 
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P(e|sJ  =  1-1  dx  dy  pvi;(x|e)  pVA'(y|J), 


(3.2-3a) 


-r-&C  -X 


.  A$-tt  ,7r+A<fc 


1  -  P(e|  6 )  = 


2  dx  dy  p(1)(x|B)  p(1)(y|l) 
-"-Ac  -X 
sL:+~  -A  4>+t 

+ J  dx  I  dy  p^(x|B)  p^(y|B) 

Ac  -r  max(-X,A<*-r) 


2  9  +  P2. 


(3.2-3b) 


(3.2-3c) 


Recall  that 


P^(x  £_) 


=  £  f„  (*+2™;D 

n=0 

oc 

’  £  Po,nP'Jx+2r"lcN-i) 

n=0 

a 

*  ■  Z  Pl,nPC<xt2F";cT-A>- 


(3.2-4) 


where  the  modulo  2r  pdf  p  ( x)  is  nonzero  for  Ix-Acj  <  it,  and  is  given  by 
(2.1-55).  For  notational  convenience,  let 


qQ(x),  c  =  cN  (unjammed) 

■ , :  »_$)  = 


(3.2-5) 


3 


With  this  notation,  we  have 


i 


iC-T  , r+i 


-d-A:  -X 


j  dy  [  (l-  ')P0jl  q0(x+2^)  +'^P1>1  qj(x+2^)] 
-x 

[('-'IPq.o  <lo(y)  +yPl,0  ql(y)J 


fL$+T 

I  dy  [(l-)2  p0s0  p0#1  q0(x)  q0(y) 


+  '  (1-^P0,0  pl,l  ql(x)  q0(y) 


+  '!  {1-)]  pl,0  p0 ,1  q0(x)  ql(y) 


+  ''2pi,i  Pl,0  qi<x)  qi(y)J 


(1"’^  P0,0P0,1  Pl(°’0)  +  YU-Y)  P0f0  Pi,i  P!(1'0) 


,^1"^  pl,0  p0,l  pi^0,1^  +  T  P1  ,lpl ,0  (3.2-6) 


From  (2.1-56),  we  recall  that 


rB 

f  dy  q^ (y)  =  F^B)  -  Fi (A)  +  \  [sgn  (B-A<)  -  sgn(A-A*)],  (3.2.7a) 

-'A 


using  the  notation 


(F0( J,  c  =  dn  (unjammed) 

F( • ;c  ,f)  = 


Fj(J,  c  =  CT  (jammed), 


.  ( 3 . 2-7  b ) 
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and  realizing  that 


Fi (x  ±  2k h)  =  Fi(x). 


Substituting  the  defined  functions  results  in 


P^i  ,j|£)  =  Fj(ii  +r  )  [Fi(A?  +  ^)  -  F^tt-A^)] 


dx  q.(x)  Fj(-x) 


+  u(i:  -  ~/2)  [Fj(A‘  +  T)  +  F.flc+ir)  -  F.(2--AO], 


in  which  the  last  term  is  zero  for  A;  <  v/2. 
The  term  may  be  written 


(3.2-7c) 


(3.2-8) 


P2  =  t(1'!)p0,0  +  '  Pl,0^ 


dy  C(l-Y)  P0>0  q0(x)  +  yPUQ  q^x)] 


>  [<1->)poi0  q0(y)  +  YP1,0  ql(y)] 


C(l-v)P0,0  +  YP1,0^  p0 ,0  P2^0,0^ 


-•  (l-^)p0j0p1)0[p2(0,l)  +  P2(1,0)]  -Y2pJ>0P2(l,l), 


(3.2-9) 
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The  probabilities  P2( 1 )  are  found  to  be 


.t-AC 


P2(i,jli)  =J  dx  qi(x)Fj(-x) 

At-r 


-  Fj (£.*-*)  [Fi -Act)  -  Fi(A<t  >t)3 


+  u(  tt/2- Act)  [F^  ( -Act)  -  (Act  -TT  )  .  F  j  (At  -n  )  ] 


where  the  last  term  is  zero  for  Ac  >  -n/2. 

The  total  P(e)  then  can  be  expressed  as 


1  1 


Pl.il  •  £  £  d-,)2-(i+i)  -,i+i  PUil.Jii) 

i  =  0  j  =0 


where 


P(eo,j  0  =  1-  P,_0  Pj>0  -  2P(>1  PJ>0  Pjli.jii) 


+  p1,0  pj  ,0  ¥'"> 


For  computation,  it  is  somewhat  more  efficient  to  use  the  following 
formulation: 


1  i 


P(e!i)  =  £  J(4j)  U-r)2-(,tJ>  V1+j  P3O.0IB) 
i =0  j=0 


(3.2-10) 


(3.2-lla) 


(3.2-llb) 


( 3 . 2- 1 2a 
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.  5  - 


r 


v- 


where 


P3(i,j:e)  -  1  -  P,,o  Pj,o  -  P,,!  Pj.o  -  Pi,0Pj,lPl(j'1li) 


+  ,0  P j ,0  ^ 


(3.2-12b) 


3.2.2  BER  Results  for  L  =  2 

The  L  =  2  bit  error  probability  for  FH/CPFSK  in  partial-band  noise 
jamming  was  calculated  as  the  average  over  bit  patterns: 

P(e)  =  ^  [ P ( e ! 1 1 1 )  +  2P ( e 1 01 1 )  +  P(e|010)],  (3.2-13) 

where  the  pattern-dependent  error  probability,  parametric  in  y,  the  jamming 
fraction,  is  given  by  (3.2-12).  The  program  listed  in  Appendix  C  was  used. 

The  results  in  Figure  3.2-2  demonstrate  that  the  BER  for  L  =  2  is 
maximized  for  a  particular  value  of  y  which  depends  on  Eb/Nj.  For  the  case 
of  Eb/NQ  =  15  dB,  h  =  0.7,  and  D  =  1.0  as  suown  in  the  figure,  we  can  observe 
that  the  worst-case  ■  value  is  approximately 


'WC 


4/(Eb/Nj)  ,  Eb/Nj  >  4  =  6  dB; 

1  »  Eb/Nj  <  6  dB; 


(3.2-14) 


for  high  Eb/Ng.  From  Figure  3.2-2  we  can  develop  the  plots  of  BER  vs. 

Eb/Nj  shown  in  Figure  3.2-3,  for  fullband  jamming  (y  =  1)  and  for  worst-case 
partial-band  jamming  (>  =  y^}.  It  is  clear  that  the  jamming  is  significantly 


more  effective  using 


WC 


than  using  y  =  1,  when  Eb/Nj  >  6  dB.  For  example, 
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a  10"4  BER  requires  E^/Nj  =  20.5  dB  when  y  =  1,  but  E^/Nj  -  37.6  dB  when 
'i  =  >w^,  a  17.1  dB  difference.  We  note  that  the  worst-case  BER  tends  to 
exhibit  an  inverse-linear  dependence  upon  E^/Nj  when  E^/Nj  >  6  dB;  that  is, 


PUC  =  -5S/(Eb/Nj).  Eb/Nj  >  6  dB, 


(3.2-15) 


for  high  E^/Nq.  This  BER  is  averaged  over  the  data  patterns;  it  is  interesting 
to  note  that  the  worst-case  BER  for  each  pattern  also  tends  to  have  the  same 
kind  of  dependence  upon  E^/Nj,  with  the  coefficients  being  approximately 
.74  for  111,  .58  for  Oil,  and  .42  for  010.  Thus  an  analysis  based  on  only 
the  "middle"  case  of  Oil  would  have  actually  represented  the  average  in  this 
instance. 

Recalling  that  Pwc  =  .21/(Eb/Nj)  for  L  =  1,  we  observe  that  the  linear 
diversity  combining  of  two  differential  phase  samples  does  not  improve  the 
system's  performance,  but  rather  degrades  it  by  about  4.4  dB.  This  difference 
is  partly  understood  from  the  fact  that,  under  a  bit  energy  constraint, 

N  =  \  ^b^O^  and  T  =  \  ^or  L  =  2.  But  this  fact  only  accounts  for  a 

3  dB  difference.  The  entire  4.4  dB  difference  may  be  attributed  to  noncoherent 
combining  losses,  which  evidently  are  so  great  for  linear  combining  of  CPFSK 
differential  phase  samples  that  even  under  an  equal  power  constraint 
(' N  =  Et/N0’  etc‘)’  E  =  2  gives  a  worse  result  than  L  =  1. 

For  example,  with  no  jamming,  E^/Nq  =  15  dB  yields  a  BER  of  1.83  *  10’10 
for  L  =  1,  and  8.3  >-  10’^  for  L  =  2;  this  difference  in  performance  is  equivalent 
to  about  3.2  dB  in  E^/hg.  So,  even  without  jamming,  the  L  =  2  combining  losses 
are  great  enough  to  give  a  worse  BER  result  for  an  equal  power  constraint. 
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3.2.3  Simplified  Calculations 

As  in  Section  3.1.3  for  L  =  1,  we  now  summarize  L  =  2  BER  results  using 
the  "typical  pattern"  approach,  rather  than  the  exact  one  involving  averaging 
over  patterns.  The  simplifying  case  is  that  of  a <p  ~  v/2,  r  =  \  =0,  and 
V  =  0,  which  gives 

r  272 

P(e;i,j)  =  1  -  pi#o  Pj,0  +  2  pilpj,0  j  dx  q i ^ x ^ F j ^ " x ^ 


tt/2 


r  72 

+  Pi,0  Pj,0  /  «*  q^xJFji-x). 

-72 


( 3 . 2- 16a ) 


where 


72 


1 


-U. (1-cosy  sinx) 

(x)  =  ^7  J  dy  cos  y  (1+Ui  +  Ui  cosy  sinx)e 
-•72 


( 3 . 2- 16b ) 


and 


r/2 


Fj(-x)  = 


dy 


■r/2 


■UJl+cosy  sinx) 

J 

1  +  cosy  sinx 


(3.2-16c) 


Using  these  expressions  with  IK  =  cp/L  =  0.405  c  ^  and  the  click 


v;; 

probabi 1 i ty 
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for  high 
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whereas  for 

apparent. 
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WC  Eb/Nj 


(3.2-17) 
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3.3  BER  COMPUTATIONS  FOR  L  >  2 

The  intricacy  of  the  FH/CPFSK  BER  calculation  was  seen  in  the  last 
subsection  to  be  considerable  for  L  =  2,  with  corresponding  amounts  of  computer 
time  needed.  In  the  following  material  we  summarize  the  extension  to  L  >  2 
of  the  direct,  or  convolutional,  analysis  approach  used  previously.  We  also 
introduce  another  numerical  approach  which  is  designed  to  control  the  amount 
of  computation  required  to  evaluate  the  BER.  Finally,  numerical  results  for 
worst-case  BER  when  L  >  2  are  presented  which  indicate  the  failure  of  linear 
diversity  combining  to  achieve  an  improvement  in  the  system  performance 
against  worst-case  partial-band  jamming. 

3.3.1  Methodology  Using  Direct  Approach 

Assuming  that  the  pdf  for  the  sum  of  L-l  diversity  samples  can  be 
calculated,  given  by 


p(L'1)(x  =  £  fJ1L"1)(x+2rn;S) 


(3.3-1) 


n=0 

in  the  notation  of  Section  3.1,  we  can  formulate  the  BER  for  the  sum  of  L 
samples  in  terms  of  this  pdf  and  the  L  =  1  pdf.  It  can  be  shown  that 


nmax  n 


P(e|£)  =  1  -  y  X 

fco  feo  m,n"m 


(3.3-2) 


where  nmax  is  the  number  of  significant  clicks,  given  by  (2.2-22b),  and 
Bm"n-m  is  the  probability  that  the  sum  is  positive  when  there  are  m  clicks 
in  the  Lth  sample  and  n-m  clicks  in  the  sum  of  the  first  L-l  samples. 


100 


■v  -v 


_  «v" 

I 


si 


■v 


3" 

9 

ft  ' 

I 


ra 


J.  S.  LEE  ASSOCIATES,  INC. 


That  is. 


,(L-1)^4>  +ir)-2Tr(n-m) 

BiL,L  -  I  *  d*  fi1)(l‘+2’m)  /  d*  ^n-m1)Cy+Z,r<n-m)3 

max[Ac--r-2irm,-(L-l)(A<t+ir)+2Tr(n-m)]  max[(L-l)(A4>  -ir)-2ir(n-m),-x] 


,(L-1)(a*+*) 


/A5+1T  / 

dX  J  d* 

max[Aj-r,2Ttn-(L-l)(A(t+7r)  ]  max[  (L-l)  (A4>-n)  ,-x+2Tin] 


For  example,  if  L  =  3  and  tt  >  t/3. 


P(ejJ)  =  1  -  Bqo  -  (B01+B1q)  -  (b02+B11+B20^’ 


where 


A:+t  r  2(A:+  h) 

dx  f^(x)  I  dy  fj2\y) 


max[2(Ao-n) ,-x] 


B0i  =  [  dx  f(1}(x)  f2(  dy  f(j2)(y) 


-x+2r 


2(A;+tv) 


■fdxf  <»(*)/  dytf’ly) 

J  -X+2tt 


.A:+r  .  2 (Ac  +n ) 

dx  f<1}(x) J  < 


dy  f^2)(y) 


2p-2a; 


=  f  dx  f!1}(x)  f 


2(Ai  ) 


dy  fi2^(y) 


(3-3-3) 


(3.3-4) 


(3.3-5a) 


(3.3-5b) 


(3.3- 5c ) 


(3.3-5d) 


( 3 . 3- 5e ) 
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’20 


/•A<f>+n 

-I  dx 

2j t-2a<(. 


/ 


2( A4+u  ) 

dy 


42)(y). 


47r-x 


(3.3-5f ) 


Notice  in  (3.3-3)  and  (3.3-5)  that  the  integration  limits  do  not  depend  upon 
m;  the  significance  of  this  fact  is  that  for  no  jatiming,  when 


P(/)(x|PN.i) 


(3.3-6a) 


and 


f(L-l) 

n 


(y)  =  p 


(L-l) 

0,n-m 


p  1)  (y(pN.^_) » 


(3.3-6b) 


then  the  click  probabilities  factor  out,  leaving  identical  integrals.  That 
is,  for  this  case 


b(D 

m,n-m 


,  DU)  n  (L'D  C 
p0,m  p0,n-m  Ln* 


( 3 . 3-7  a ) 


where 


<Av+  • 


Cn 


dx  pU)(xk'N,i) 


,L-1)  (A4>-Ht) 

dy  P^L"1}(y|pN*A) 


max[Ac-r  ,2m- (L-l)  (a^+tt)  ]  max  [(L-l)  (a<j>-tt)  ,2irn-x] 


numerically,  the  integrals  have  to  be  computed  once,  rather  than  n  +  1  times. 
However,  for  the  general  case,  the  click  probabilities  do  not  easily  factor 
out  of  f^  and  fj^”^  (see  (3-2-4),  for  example],  and  the  double  integral 
must  be  computed  for  each  m.  In  view  of  the  increasingly  complex  and  time- 
consuming  computations  required  by  the  direct  method  as  the  order  of  diversity 
L  increases,  we  have  developed  numerical  methods  for  implementing  a  character¬ 
istic  function  or  transform  approach. 
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3.3.2  Methodology  Using  Transform  Approach 

As  we  have  noted,  it  is  sufficient  for  BER  calculations  to  model  the 
single-sample  pdf  as  having  a  finite  maximum  number  of  clicks: 


nmax 


P{1)(xjs)  =  fn(x+2Tin;6). 


(3.3-8) 


n=0 


This  model  is  sufficient  because  numbers  of  clicks  greater  than  nmax  simply 
do  not  necessarily  enter  into  the  BER  calculation.  Given  the  pdf  (3.3-8),  the 
range  of  the  differential  phase  (x)  which  is  significant  then  is 


:-r-2"nmax  <  x  <  M+r., 


(3.3-9) 


Since  the  pdf  is  zero  outside  this  range,  we  can  consider  numerical  evaluation 
of  the  characteristic  function,  by  means  of  the  DFT  (discrete  Fourier  transform) 
of  the  pdf.  Theoretically,  it  is  well  understood  that  if 


p 

r  a 

v'(v)  -  j  dx  e-'-’2'*  p^(x)  =  # 

P(1)(x)| 

(3.3-10) 

1 

*y. 

—  OC 

ft 

is  the  characteristic  function  for  one  sample,  then 

I  I  L  I 

*(v)|  = 


P^(x) 


(3.3-11) 


that  is,  the  Lth  power  of  <  is  the  characteristic  function  for  the  sum  of 
L  independent  and  identically-distributed  samples.  The  challenge  is  to  utilize 
DFT  methods  and  parameters  which  will  provide  sufficient  accuracy  to  evaluate 
the  BER. 
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3. 3. 2.1  DFT  Size  Considerations 

Figure  3.3-1  illustrates  the  fact  that,  since  the  L  =  1  pdf  for  a  given 
pattern  is  nonzero  for  the  range  A^-^immax  <  x  <  A<j>+ir,  the  L-hop  pdf  is 
nonzero  for  the  range  L(a<±  - tt-2  nmax)  <  x  <  L(A<j>+n).  Therefore,  to  avoid 
aliasing,  the  N-point  DFT  must  consider  pdf  samples  on  an  interval  at  least 
as  long  as 

X  =  Nax  >  L(nmax+1)27:  radians.  (3.3-12) 

From  (2.2-22b),  a  conservative  estimate  of  nmax,  the  number  of  "significant" 

2 

clicks,  is  nmax  =  L-l.  Thus  the  interval  should  be  at  least  L  2 t  radians  in 
length.  At  the  same  time,  it  is  convenient  to  stipulate  that  in  ~  radians 
there  are  exactly  an  integer  number  of  sampling  intervals,  that  is, 

~  =  integer  =  ng.  (3.3-13) 

This  implies  that 

X  =  Nax  =  N  •  T/nQ  >  2*l2;  (3.3-14a) 

or  that  the  DFT  size  N  is  bounded  by 

N  >  2rL2/Ax  =  2L2nQ  .  (3.3-14b) 

For  example,  if  we  specify  that  ng  =  128  samples  are  desired  over  *  radians 
(ax  =  ~/128),  then  the  following  DFT  sizes  are  required: 


Vi  r  r  *  *  4"  c 

>1^5  %  S  ^  >  VV 
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minimum  N 


referred  N 


In  this  table,  we  choose  a  "preferred  N"  value  that  is  the  next  higher  power 
of  two,  in  order  to  provide  "zero  fill"  for  interpolating  the  characteristic 
function  with  more  closely-spaced  points  in  the  transform  domain.  It  seems 
unlikely,  however,  that  there  is  any  "best"  window  or  window  length  in  this 
application  of  the  DFT,  since  the  effective  width  (standard  deviation)  of  the 
pdf  will  vary  with  SNR,  and  it  is  impractical  to  try  to  match  the  two  widths 
in  some  sense  for  each  SNR  value. 

3. 3. 2. 2  Formulation  of  BER  Using  DFT 

The  strategy  we  have  adopted  for  calculating  the  BER  is  designed  to 
exploit  the  fact  that  a  finite  number  of  clicks  is  involved  in  calculation 
of  the  probability  of  a  correct  decision,  P(C)  =  l-P(e).  Thus  our  approach 
is  to  utilize  DFT  methods  to  calculate 


P(C)  =  j  [P(C  ! Ill )  +  2P(C  |011)  +  P(C  |010)3, 


(3.3-15) 


where  the  conditional  probability  of  a  correct  decision  is 

,l(ii+") 

P(Cji)  -  J  dx  p^(x|j?) 

0 

,jL{L:* dx  j”  dv  P2"*  ’ 

0  -co  L  A<f-r-2Trnmax 


(3.3-16) 


lira 


we  use  Simpson's  rule  to  perform  the  outer  integral,  with  the  L-diversity 
pdf  approximated  by 

-N-l  _  N-l  - 

p(L>(m/.x)  -  av  Y.  ej27"rixiv[iy  Y.  e'J2rrntySvP(niJ')|  L  ,  (3.3-17) 


in  which  p  is  the  periodic  extension  of  pv  1 .  Since  x  =  Ay  =  X/N  and 
L\  =  1/X,  then  :  x  •  Ay  =  1/N  and  (3.3-17)  becomes 


p(L)(mix)  ~  XL_1  IDFT  [DFT(p)]L  . 


(3.3-18) 


The  use  of  the  inverse  DFT  (IDFT)  in  (3.3-17)  and  (3.3-18)  is  based 


on  the  fact  that  the  finite  Fourier  transform  of  the  periodic  extension  of  the 
(finite  domain)  L  =  1  pdf  is  (a)  equal  to  the  pdf's  characteristic  function 
at  the  points  =  n  •  and  (b)  periodic  with  period  1/ax  when  the  finite 
Fourier  transform  is  approximated  by  a  DFT.  That  is, 


=  P  (x) 


dx  e"'*2,  '*  ^ 


p^  (x+kX) 


f  j2r'x  p(1)(x)  +  ej2rvX  f°dx  e-j2"vxp(1)(x) 

-A/ 


^vl'(  O.  v  =  n  •  Av  =  n/X; 


(3.3-19) 
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and  if 


%  (v)  »  DFT  {  p(x)  }  , 


rN  C(n+N)Av]  = 


(3.3-20a) 


bo  : 


( 3. 3-20b) 


Therefore  approximating  the  inverse  Fourier  transform  integral  by  a  summation 
from  .  =  -  N :./2  to  .  =  Niv/2  is  the  same  as  a  summation  from  v  =  0  to 
'■  =  ( N- 1 )  -  . 

The  approximation,  besides  utilizing  a  discrete  sum  to  calculate 
an  integral,  is  predicated  on  the  characteristic  function's  vanishing  for 
-  N../2;  otherwise,  it  will  be  distorted  due  to  aliasing.  This 
vanishing  never  exactly  occurs,  since  the  pdf  is  nonzero  over  a  finite 
interval.  However,  the  more  smoothly  the  pdf  decreases  at  its  end  points, 
the  less  "bandwidth"  is  required,  and  aliasing  of  the  characteristic 
function  is  minimized.  For  example,  the  least  "smooth"  pdf  occurs  for  zero 
SNR,  and  is  given  by 


nmax 

p^(x;.  =0)  =  pn[u(A;  -r+Znr  )  -  u(^ -hr+2n  *)]  , 


(3.3-21) 


in  which  u(*)  is  the  unit  step  function  and  the  {  pn  }  are  click  probabilities. 
The  magnitude  of  the  corresponding  characteristic  function  is  bounded  by 
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for  v  =  Nav/2  and  A x  =  t/128. 


(3.3-22) 


i 
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This  is  26  dB  down  from  the  maximum  value  of  v>(0)  =  1.  A  triangular  pdf 
for  no  clicks  would  be  52  dB  down  at  the  N/2  folding  point,  and  a  Gaussian 
density  with  o-  =  1/:  would  be  down  the  following  amounts  for  different  SNR 
values,  c  ,  and  :  x  =  it/128: 


(expression) 
0  dB 
10  dB 
20  dB 


Attenuation 
Bound  at  NAv/2 


0.54/; (ax)2  (dB) 
901  dB 
90  dB 
9  dB. 


3. 3. 2. 3  Normal izations 

Referring  now  to  equations  (3.3-16)  and  (3.3-17),  the  characteristic 
function  for  L  =  1  is  approximated  by 


v'U)(„  )  = 


J.  +7 


3-j2..vx  p(l)(x)  r  x  OFT  { p }. 


dx  e 

L: -r-?7nmax 

As  a  means  for  controlling  the  approximation  error,  we  utilize  the  fact 


(3. 


3.23) 
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niiidx 

■  I  [(l->)p0,n+vp1>n]. 


(3.3-24) 


Therefore  for  each  intersymbol  .nterference  pattern  we  adjust  the  DFT 
values  by  a  normalization  factor  to  produce 


DFT(k)'  ■  ^  , 


(3.3-25a) 


where 


NF  =  X  •  DFT(0)/V'  1  (0) 


(3.3-25b) 


3. 3. 2. 4  Averaging  Over  Data  Patterns 


To  minimize  the  number  of  inverse  DFT  operations,  we  perform  averaging 
over  data  patterns  prior  to  the  inverse  DFT.  That  is,  the  pdf  for  the  diversity 
sum  is  computed  as 


p^ L  ^  (nr  x) 


=  X1'1  IDFT  |  [DFT  (p,lll)/NF(lll)]L 
+  \  [DFT  (p,011)/NF(011)]L 
+  l  [DFT  (p,010)/NF(010)JL 


(3.3-26) 


3.3.3  Results  Using  DFT  Method 


Since  we  previously  have  calculated  the  FH/CPFSK  bit  error  probability 


for  L  =  2  by  the  direct  method,  the  accuracy  of  the  DFT  method's  results 
may  be  discerned  by  making  comparisons  for  this  case.  Also,  the  values  of 
the  normalization  factors  for  the  different  patterns  give  a  general  indication 
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of  how  well  the  discretized  pdf  is  representing  the  actual,  continuous 
distribution. 

The  worst-case  BER  was  found  by  computing  a  few  points  on  BER  vs  y  curves 
like  the  ones  shown  in  Figure  3.2-2,  just  enough  points  to  establish  the  P(e) 
maximum  with  respect  to  > .  The  program  included  in  Appendix  D  was  used,  and 
each  point  required  between  five  and  ten  minutes  of  computation,  depending 
on  L.  Normalization  factors  typically  were  between  0.95  and  1.05,  indicating 
that  numerical  integration  of  the  L  =  1  pdf  using  2r/ox  =  256  points  and  a 
rectangular  rule  would  yield  about  a  5%  error,  if  the  factors  were  not 
employed. 

3. 3. 3.1  Results  for  Eb/NQ  =  15  dB 


Figure  3.3-2  summarizes  the  worst-case  partial-band  noise  jamming 
performance  of  FH/CPFSK  using  discriminator  detection  and  linear  diversity 
combining.  The  curves  were  drawn  using  the  data  in  Table  3.3-1.  We  have 
already  noted  in  Section  3.2  that  the  L  =  2  performance  is  about  4.4  dB 
worst  than  for  L  =  1  (no  diversity);  these  previous  results  of  direct 
calculation  are  included  in  the  figure  for  reference.  Worst-case  BER  results 
for  1=2  using  the  DFT  method  are  also  presented  for  the  portion  of  the  graph 
lying  between  10  and  30  dB.  Pictorially,  the  L  =  2  results  using  the  two 
different  methods  are  barely  distinguishable;  the  data  in  Table  3.3-1  indicates 
that  the  DFT  method  gives  a  BER  about  3%  high,  relative  to  the  direct  method. 
This  was  considered  to  be  an  acceptable  degree  of  agreement  between  the  two 
methods,  so  that  further  development  of  the  computer  program,  such  as  more 
elaborate  normalizations,  was  not  pursued. 
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The  L  =  3  results  in  Figure  3.3-2  reflect  a  2.4  dB  worse  performance 
than  for  L  =  2,  when  the  jammer  power  is  strong  (low  E^/Ng).  As  E^/Nj 
increases,  the  BER  for  this  number  of  hops/bit  converges  to  an  asymptote  of 
about  4  *  10'4,  the  performance  for  L  =  3  and  no  jamming.  Evidently  the 
unjammed  noncoherent  combining  loss  at  E^/Ng  =  15  dB  is  about  5.2  dB,  because 
(from  Figure  3.1-3)  a  4  *  10~4  BER  for  L  =  1  occurs  for  E^/Ng  :  9.8  dB.  This 
loss  contrasts  with  the  amount  in  strong  jamming,  about  6.8  dB. 

As  anticipated,  the  L  =  4  results  in  Figure  3.3-2  follow  the  trend  of  the 
BER  increasing  with  L,  at  least  for  strong  jamming.  We  observe  that 
for  15  dB  <  E^/Nj  -  25  dB,  the  L  =  4  worst-case  performance  dips  below  that  of 
L  =  3,  before  settling  to  its  unjammed  value  of  9  ^  10"4  (5.7  dB  unjammed 
noncoherent  combining  loss).  This  phenomenon  is  very  interesting,  since  we 
are  seeking  a  BER  behavior  which  decreases  with  L.  But  in  this  figure,  we 
cannot  discern  the  trend  because  the  L  =  3  and  L  =  4  performances  are  influenced 
so  much  by  thermal  noise. 

3.3  Results  for  E^/Ng  =  20  dB 

n  figure  3.3-3  we  show  L  =  2,  3,  and  4  BER  results  for  E^/Ng  =  20  dB, 

plotted  from  the  data  in  Table  3.3-2.  For  each  diversity  value  (L),  the 

thermal  noise  is  not  influential  -  that  is,  the  unjammed  BER  is  much  less  than 
-5 

10  -  so  we  can  observe  the  trend  of  the  relative  performances  for  L  =  3 

and  L  =  4.  What  we  see  is  the  fact  that,  for  negligible  thermal  noise,  the 
worst-case  BER  performance  for  L  =  4  lies  between,  and  parallel  to,  those  for 
L  =  2  and  L  =  3.  Thus,  although  a  slight  improvement  is  made  (about  1  dB), 
there  is  no  "diversity  gain"  improvement  for  higher  L.  If  there  were  such  an 
improvement,  the  negative  slope  of  the  (log)  P(e)  curve  vs  E^/N j  in  dB  would 
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FIGURE  3.3-3  WORST-CASE  FH/CPFSK  BER  VS  E^Nj  FOR  E^/Nq  =  20  dB  AN 
THE  NUMBER  OF  HOPS/BIT  (L)  VARIED 
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TABLE  3.3-2  DATA  FOR  FIGURE  3.3-3 


p 

Kwc 

^WC 

p 

KWC 

Ywc 

p 

KWC 

YWC 

0.329 

1.0 

0.37 

1.0 

0.385 

1.0 

0.173 

1.0 

0.218 

1.0 

0.263 

1.0 

i 

0.17 

1.0 

5.5(2) 

0.3 

8 . 6  ( -2 ) 

0.48 

8.9(-2) 

0.55 

1 . 74 ( - 2 ) 

9.5(-2) 

2.8(-2) 

0.15 

2. 6(-2) 

0.16 

5 . 5  ( -  3 ) 

3 . 0  ( -  2 ) 

9.0(-3) 

4.8(-2) 

7 . 8 ( -3 ) 

5 . 0  ( -2 ) 

2.8(-3) 

1.5(-2) 

5 . 6( -4 ) 

3 . 0( -3 ) 

9 . 0  ( -4 ) 

4.8(  -3) 

7.8(-4) 

5 . 0  ( -  3 ) 

5. 5( -5) 

3 . 0  ( -4 ) 

9 . 0  ( -  5 ) 

4.8( -4) 

7 . 7  ( -  5 ) 

5.0(-4) 

5.5(-6) 

3 . 0  ( -  5 ) 

9 . 0 ( -6 ) 

4 .8  ( -  5) 

7 .7 ( -6) 

5 . 0  ( -  5 ) 
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be  greater  than  unity. 

3. 3. 3. 3  Explanation  of  the  Observed  Diversity  Behavior 

What  then  accounts  for  the  slight  worst-case  BER  improvement  going  from 
L  =  3  to  L  =  4?  The  apparent  explanation  can  be  made  with  the  help  of  Figures 
3.3-4  to  3.3-7.  In  Figure  3.3-4,  we  show  the  pattern-averaged  pdf's  for  a 
single  hop  sample  and  for  the  sum  of  two  hop  samples  when  L  =  2  for  particular 
values  of  "i  ,  E^/Nq,  and  E^/Nj.  The  linear  scale  is  such  that  the  click 
contributions  are  not  noticeable,  so  in  Figure  3.3-5  we  repeat  the  same 
information,  using  a  logarithmic  scale.  The  superposition  of  the  three  pattern- 
dependent  pdf's  constituting  the  averaged  pdf  for  L  =  1  is  now  evident  in 
Figure  3.3-5,  and  we  call  the  reader's  attention  to  two  properties  of  the 
two-sample  pdf:  (1)  its  major  peak  is  shifted  to  the  right,  located  at  &C  , 
compared  to  L:  for  one  sample;  (2)  the  width  of  its  lobes  are  increased  over 
that  for  one  sample.  Thus,  while  for  L  =  2  the  pdf  shifts  to  the  right,  tending 
to  decrease  the  error,  it  also  is  spreading  out,  and  furthermore  the  "mass"  under 
the  one-click  lobe  is  greater  than  for  one  sample;  both  these  latter  trends 
tend  to  increase  the  error. 

In  Figures  3.3-6  and  3.3-7  we  present  similar  pdf  illustrations,  for 
L  =  3  and  L  =  4,  respectively.  The  trend  is  for  the  nonzero  click  lobes  to 
become  more  significant  as  L  increases,  since  the  SNR  is  decreasing  for 
constant  bit  energy  (c^  =  ( E^/Nq ) /L ) .  However,  as  L  goes  from  3  to  4,  note 
that  the  peak  of  the  1-click  lobe  crosses  zero.  This  means  that  for  high  SNR 
and  L  =  4,  a  correct  decision  is  made  even  when  there  is  a  single  click, 
whereas  for  L  =  3  a  correct  decision  is  not  made  if  there  are  any  clicks  at  all. 


PATTERN-AVERAGED  DIFFERENTIAL  PHASE  SUM  PDF'S  FOR  L  =  2  (LINEAR  SCALE) 
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Thus  L  =  4  has  a  slight  advantage  over  L  =  3,  but  does  not  offer  a  "diversity 
gain"  improvement. 

3. 3. 3. 4  Significance  of  the  Diversity  Behavior 

The  BER  results  shown  in  Figures  3.3-2  and  3.3-3  demonstrate  the  fact 
that  a  diversity  gain  against  worst-case  partial-band  noise  jamming  is  not 
realized  for  the  linear  combining  of  FH/CPFSK  differential  phase  samples  from 
different  hops.  The  same  fact  is  true  for  FH/BFSK,  but  it  was  plausible  to 
conjecture  that  the  nonlinear  demodulation  procedures  for  CPFSK  (e.g.,  the 
bandpass  limiter)  would  provide  the  sort  of  nonlinear  processing  that  produces 
a  diversity  gain.  Now  that  this  conjecture  has  been  disproved,  we  observe  with 
hindsight  that  the  limiter  (ideally)  has  no  effect  on  the  phase  of  the  signal, 
and  that  there  is  no  mechanism  in  the  limiter/discriminator  receiver  which  acts 
to  de-emphasize  samples  from  jammed  hops. 

It  is  easy  to  show  that  some  form  of  diversity  gain  is  possible  with  the 
appropriate  processing,  at  least  for  low  thermal  noise.  Consider,  for  example, 
the  scheme  in  which  it  is  assumed  that  perfect  side  information  is  available  on 
which  hops  are  jammed,  and  samples  from  jammed  hops  are  excluded  from  the 
diversity  sum  unless  all  hops  are  jammed.  Then,  for  no  thermal  noise,  no 
error  occurs  unless  all  hops  are  jammed: 


I 

urn, 


& 

$ 


r'J 

XN 


-v- 


s. 


P(e)  =  \  L  gL(>  Rj),  Rj  £  Efa/Nj, 


(3.3-27) 


where  gL(  •)  is  the  BER  vs.  E^/Nq  function  for  the  sum  of  L  samples.  Differen¬ 
tiation  with  respect  to  yields  the  equation 
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L  gL(x)  +  x  gL(x)  =0,  x  =  yRj. 


(3.3-28) 


Clearly,  if  this  equation  is  solved  for  x  =  x^,  then  the  worst-case  partial- 


band  jamming  fraction  is 


Rj  •  RJ  >  V 


1 »  RJ  <  XL" 


(3.3-29) 


Substituting  the  value  of  ^ur  back  into  (3.3-27),  we  find  that 


x  \*" 

—  )  9L(\),  Rj  >  \\ 
RJ  / 


gL(Rj) *  Rj  < 


(3.3-31) 


That  is,  Pwc  is  proportional  to  RjL  for  Rj  >  xL.  On  a  log  P(e)  vs  Eb/Nj  ( dB ) 
plot,  Pwc  is  a  straight  line  with  slope  -L,  tangent  to  the  Y  =  1  error  curve 
at  Rj  =  x^. 

Using  data  from  Tables  3.3-1  and  3.3-2  for  y  =  1,  and  constructing  tangents 
with  slopes  equal  to  -L,  we  obtain  the  predicted  ideal  diversity  performances 
for  FH/CPFSK  in  worst-case  PBNJ  as  shown  in  Figure  3.3-8.  The  higher  L  curves 
must  eventually  cross  those  for  lower  L  because  of  their  greater  slope,  and 
indeed  do  as  shown.  The  figure  illustrates  that  a  diversity  gain  is  realized 
for  this  ideal  situation,  and  if  the  optimum  value  of  L  is  always  used,  the 
worst-case  BER  can  be  made  to  approach  within  5  or  6  dB  of  the  unjammed  CPFSK. 
performance,  compared  to  over  30  dB  without  diversity. 


FIGURE  3. 
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In  practice,  thermal  noise  is  always  present,  giving  rise  to  noncoherent 


combining  losses  which  prevent  achievement  of  the  ideal  diversity  gains 


pictured  in  Figure  3.3-8.  However,  these  results  do  suggest  that  some  form 


of  combining  the  differential  phase  samples  from  the  different  hops  will 


accomplish  a  diversity  improvement  of  some  degree. 
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4.0  EVALUATION  OF  DIVERSITY  PERFORMANCE  USING  DIFFERENTIAL  DETECTION 

An  alternate  receiver  configuration  for  FH/CPFSK  is  one  in  which  the 
limiter-discriminator  is  replaced  with  a  differential  detector.  In  this 
section,  we  find  the  BER  obtained  by  FH/CPFSK  in  partial-band  noise  jamming 
when  a  differential  detector  is  used. 

4.1  ANALYSIS  FOR  NO  DIVERSITY 

The  analysis  of  differential  detection  of  narrowband  FM  presented  by 
Simon  and  Wang  [11]  relies  on  the  phase  distribution  theory  of  Pawula,  Rice, 
and  Roberts  [8],  applied  to  the  binary  FM  communications  problem  by  Pawula  [3], 
The  derivation  of  the  error  probability  using  this  theory  is  somewhat  involved. 
In  this  section,  we  present  a  simple  derivation  of  the  binary  FM  bit 
error  rate  (BER)  using  differential  detection. 

The  differential  detector,  shown  in  Figure  4.1-1,  develops  the  output 


z(t)  =  |  R(t)R(t-T)sin[e(t)-6(t-T)] 


(4.1-1) 


from  the  narrowband  waveform 


x(t)  =  R(  t)cos[-'ot  +  e(t)] 


=  x  (t)cos^ot  -  x  (t)sinwot. 


(4. l-2a) 


(4. 1 - 2b ) 
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The  rationale  for  using  this  detector  is  that,  for  binary  FM  or  CPFSK 
(continuous  phase  frequency-shift-keying),  the  bit  information  is  contained 
in  the  sign  of  the  difference  between  the  signal  phase  values  at  bit  sampling 
times  t  =  kT.  The  sign  of  the  q  antity  z  at  the  sampling  times  therefore  can 
be  used  for  the  bit  decision.  Without  noise  x(t)  is  the  signal 

s(t)  =  a(t)/25  cos[u)0t  +  4>(t)]  ,  f 4. 

where  S  is  the  signal  power,  a(t)  is  amplitude  modulation  induced  by  the 

receiver  filtering,  and  c(t)  is  the  information  phase  modulation  waveform 

after  distortion  due  to  receiver  filtering.  Prior  to  this  filtering,  the 

signal  has  constant  amplitude  and  the  bit  information  is  conveyed  by  the 

instantaneous  frequency,  which  is  either  fn  +  f  .  for  a  "mark"  or  f„  -  f  . 

u  a  u  a 

for  a  "space".  Therefore,  neglecting  distortion,  the  phase  difference  at 
=  kT  is 


=  :(tk)  -  :(tk-T)  =  ±rh. 


where  h  =  2fdT  is  the  modulation  index. 


Our  approach  to  analyzing  the  BER  for  the  differential  detector  recognizes 
that  the  random  variable  z  is  a  quadratic  form  in  the  quadrature  components 
xc(t)  and  x  (t): 


Z  =  4  [*,( t)x  ( t-T)-x  (t)x.(t-T)] . 


i 
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With  this  understanding,  as  shown  below  we  can  apply  a  convenient  equivalence 
for  random  variables  of  this  type. 

4.1.1  A  Statistical  Equivalence 

Let  the  four-dimensional  vector  of  random  variables  x.  =  (x1 


matrix 


Kx(n,C)  = 


an. 

with 

mean  vector 

— x  =  <ml 

,m2 ,mg ,m^ 

ct 

0 

no  j  o2 

Coio2 

0 

-£oi°2 

no  i  a2 

ne 

1  °2 

-E,olc2 

0 

ic 

=3 

Q 

Q 

fo 

0 

Q>\ 

(4.1-6) 


We  use  aj  to  denote  the  transpose  of  the  column  vector  a^.  In  [12]  and  [13]  it 
is  shown  that  the  quadratic  form 


y  =  2  (xlx3  +  X2X4} 


(4.1-7) 


is  equal  to  the  difference  of  two  scaled,  independent  noncentral  chi-squared 
random  variables  with  two  degrees  of  freedom,  denoted  by 

y^c1y/(2;d1)  -  c2x2(2;d2),  (4.1-8) 

where  d.^  and  d2  are  the  noncentrality  parameters.  In  terms  of  the  components 
of  mx  and  Kx,  the  parameters  are 


cl,2  =  ±  *0 


2  2  2  2  -- 

o^(m^+m2)+  o^m^+m^)  +  2oiQ2|/l"C2(n^nym2m^)  -  2o)g2^(mim4~m2m; 


( 4 . l-9a ) 

(4. l-9b) 


2o^o^  /l- 12  ( /1-S2  n) 
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4.1.2  Application  of  the  Equivalence 

Comparing  (4.1-5)  and  (4.1-7),  we  define 


? 

i* 


s 


xc(t) 

Xs(t) 

-x$(t-T) 

xc(t-T) 


(4.1-10) 


In  this  case,  the  mean  vector  m  is 


a(t)/25  cos<»(t) 
a(t)/?5  sin<j>(t) 
-a(t-T)/25  sin  4>  (t-T) 
a(t-T)/2S  cos  4  (t-T) 


Aj  cos<h 
A1  sin4>i 


-A2  si n4>2  • 

A2  C0S4>2 


>2-  b 


(4.1-11)  £ 


For  stationary  bandpass  Gaussian  noise,  of  =  of  =  NqW^,  where  WN  =  Wjp  is  the 
noise  bandwidth  of  the  receiver  filter,  and  the  autocorrelation  function  is 

Rn(x)  =  NQWN[r(T)  cosu)0t  +  X(T)sino)0f];  (4. 

thus  c  =  r(T)  e  r,  and  n  =  -a(T)  s  -a.  The  cross-quadrature  correlation 
coefficient  A  is  zero  if  the  filter  passband  is  symmetric  about  the  center 
frequency  fg. 


(4.1-12)  & 
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Substituting  appropriately  in  (4.1-8)  and  (4.1-9),  then,  we  find  that  the 


differential  detector  output  has  the  distribution 


z~c1x2{2;d1)  -  c2x2(2,d2) 


Cl,2  =5N0Mn(^T  X) 


(4.1- 13a ) 


(4.1- 13b ) 


Pi  +  P2  —  2/^Fi  ^l-r2sin  -2 
A-r2  ( A-r2  +  a) 


7p  i p  ?  r  cos  A4> 


(4.1- 13c ) 


In  ( 4 . 1- 13c )  we  define  SNR's  ^  =  A?/2NqWn  for  i  =  1,2.  Further,  defining 

A  A  --  . 

U  =  (pi+c:)/2  and  W  -  ^p1p£,  we  have 

2(U  -  r  Wcosa$  ±  /l-r^W  s i n A<t> } 

^^2  =  /  —  _  , 

A-r2  (/  1-r2  +  A) 


(4.1-14) 
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4.1.3  BER  for  L=1  Without  Jamming 

Now  since  the  quantities  U,  W,  and  &4>  are  affected  by  intersymbol  inter¬ 
ference,  the  BER  must  be  averaged  over  possible  bit  patterns.  Using  the 
symmetry  of  complementary  patterns,  the  BER  can  be  written 


3 


«Li» 

t 

,VJ 


P(e;r,x)  =  Pr{pattern  j} 

0 

x  |  [Pr  {z  <  0  |  A<j>j  ,Uj  ,Wj  } 

+  Pr {z  >  01  -A^.U  Wj)] 

=  Y.  Pr{Pattern  P(e;r,x(A4>j,Uj,Wj).  (4.1-15) 

J 

In  view  of  ( 4 . 1- 13a )  the  probability  that  z  is  less  than  zero  is  known  to  be 
expressible  in  terms  of  Marcum's  Q-function  and  the  Iq(-)  Bessel  function  [14]: 


v;* 


•N 


Pr{z  <  0|a*,U,W]  =  \  [1  -  Q(/^,/aJ)  +  Q(/aJ,/EJ)] 

«"  V 

1  1  Cl’C2  rrn  Vbl 

2  Z1+c2  exp  2 

) » 

(4.1 - 16a ) 

& 

“A* 

where 

t  C2^2^c^  1  r - 

a.  =  -  —  - -  =  - ?  {  U  -  rWcosAcf  -  /l-r2W  sinA<j>} 

1  C1  c2  1 -TC 

(4. 1- 16b) 
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• 

• 

h  - 

D1  “  cY+cz 


— ^  {  U  -  rWcosA<j>  +  v'T-r^W  sinA<|>} 
1-/ 


( 4 . 1- 16c ) 


Similarly, 


Pr(z  >  0  |-A<j>,U,W} 

=  \  [1  -  Q^,/^)  +  Q(/^,v^)] 

i  ___ 

+  2  cj+c^  eXp{  “  }  I0^a0^)0^  * 


(_  Cjdj(-A$)  ) 


0  cl+c2 


^  A  c2^2^"a<^^ 
bO  c^+Cg 


Cl+C2 


Cjd^A^) 

Ci+c2  _  bl' 


(4. l-17a) 


(4.1- 17b ) 


(4.1-17C) 


Because  a^  =  a^  and  bp  =  b^,  the  terms  containing  the  exponential  in  (4.1-16a) 
and  (4.1-17a)  cancel  when  they  are  averaged  to  get  the  conditional  BER,  resultiny  in 


P(e;  r,x |a$,U,W)  =  (1  -  Q(/b^,/a^)  +  Qf^./b^)  J 


(4.1-18) 


4. 1.3.1  Comparison  with  Other  Analyses 

Note  that  the  expression  (4.1-18)  does  not  depend  on  the  cross-quadrature 
correlation  coefficient  so  this  quantity  can  be  assumed  equal  to  zero  for 
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convenience.  This  same  kind  of  symmetry  (for  equally  probable  data  symbols) 
has  been  noted  for  DPSK  [13],  [16]  with  respect  to  r,  the  same-quadrature 
correlation  coefficient. 

The  conditional  binary  FM  error  probability  given  by  (4.1-18)  can  be 
shown  [7,  p  85]  to  be  identical  with  the  expressions  given  in  [11]: 


P(e;  r.olM,  U,  W)  »  |  [1  -  /l-B2/*2  Ie(B/a,a)], 


(4. 1 - 1 9a ) 


where  Ie(*,*)  is  the  Rice  function,  with 


a  =  Uj+bjJ/2  ,  8  =  /ajbj. 


Calculations  of  (4.1-19)  involve  numerical  integration,  since  an  equivalent 


expression  is  [11] 


jr 

P(e;  r,  0 |a$,  U,  W)  =  ^ 


f /’ 


pr  -(g-Bcose; 
a-B  cose 


(4.1-20) 


Accurate  calculation  of  the  BER  using  (4.1-18)  involving  Marcum's  Q-function 
is  considered  easier  because  of  the  possibility  of  a  singularity  in  the  inte¬ 
grand  of  (4.1-20) . 

It  is  noteworthy  that  the  error  expression  (4.1-19a)  represents  the 
general  form  for  the  BER  for  any  modulation  scheme  for  which  the  receiver  output 


at  the  sampling  instant  can  be  written  as 


z  =  R  cos(c i  ±  $2) ,  R  51  0, 


(4.1-lSb) 


(4.1-21)  ^  \ 


with  -1  and  c2  distributed  as  the  phase  of  a  sinusoid  in  noise.  This 
observation,  due  to  Jain  [17],  is  applied  by  him  (with  our  parameters  r=x=0 
and  £v=ty2)  to  undistorted  detection  of  (a)  hard-limited  PSK  with  a  perfect 
reference,  (b)  PSK  with  a  noisy  reference,  (c)  DPSK  with  differential  detection 
(pi=pJ»  (d)  binary  FM  with  discriminator  detection  and  without  integrate-and- 
dump  output  filtering  (oi?cz),  and  (e)  BFSK. 

4. 1.3.2  A  Useful  Approximation 

The  noncentral  chi-squared  distribution  with  v  degrees  of  freedom  and 
noncentrality  parameter  d  is  well  approximated  by  ([18], [19]) 

Pr ; Y  >  x\,d}  =  Q(>  x  -(v-l)/2  -  /d  +  “(v-D/2),  (4.1-22) 

where  Q(-)  (with  one  argument)  is  the  Gaussian  complementary  distribution 
function.  Therefore,  Marcum's  Q-function  is  approximated  by 

Q(-,0  =  Pr(x:  >  £2  l2,cx2 } 

=  USeTTffi  -  s* 2  +  1/2 ) 

=  QU  -  a). 

The  conditional  BER  (4.1-18)  then  is  approximated  by 
P(e;r,0  jAi.U.W)  -  Q( -  /aT  ) . 
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For  differential  detection  of  narrowband  digital  FM,  averaging  the 
conditional  BER  over  bit  patterns  gives  the  unconditional  error  probability 


P(e)  =  |  { [P(ejOOO)  +  P(e| 111)] 


+  2[P(e | Oil )  +  P( e [ 100) ] 


+  [P(ejOlO)  +  P(e 1 101 ) ]  }  . 


(4.1-25) 


To  evaluate  the  accuracy  of  the  approximation  given  by  (4.1-24),  and  to 
compare  the  exact  formula  (4.1-18)  with  the  results  presented  in  [11],  we  use 
the  sample  case  of  h  =  0.7,  D  =  1.0,  and  A  =  0.  The  exact  values  for  U,  W, 
and  1:  needed  to  calculate  a1  and  b1  according  to  (4.1-16)  were  given  in 
Table  2.1-3,  and  the  same-quadrature  correlation  coefficient  equals 
e"‘  =  .0432.  Since  in  [11]  the  parameters  for  the  Oil  and  100  patterns  are 
approximated,  we  shall  use  the  approximate  formulas 


U  (Oil)  =  i  [U( 111)  +  U(010) ] , 


W  (Oil) 


U(lll)  U(010), 


(Oil)  :  \  (111)  +  A4(010)] 


The  BER  for  this  example  case  is  shown  plotted  against  E^/Nq  in  Figure  4.1-2. 
These  results  indicate  that  the  approximation  tends  to  give  a  low  estimate  of 
the  BER,  but  is  quite  good  for  high  SNR. 


K 


ft  \ 


(4. l-26a) 


(4.1-26b) 


(4.1-26c) 


■  »>  ,  r 


.  - .  S 

«L 
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4.1.4  BER  for  L  =  1  With  Jamming 

The  case  of  partial -band  noise  jamming  for  L  =  1  hop/bit  can  be  treated 
as  an  extension  of  the  unjammed  case  in  the  following  way: 

P(e;  y)  =  (1  -  t)P3  +  fPb,  (4.1-27) 

in  which 

Pa  =  pattern-averaged  BER  for  CNR  =  (4.1-28a) 

P.  =  pattern-averaged  BER  for  CNR  =  E ^/N-j. ,  ( 4.1-28b ) 

where  the  averaging  is  of  the  conditional  BER  given  by  (4.1-18)  over  the 
pattern-dependent  quantities  U,  W,  and  a <p  as  listed  in  Table  2.1-3. 

Numerical  calculations  of  the  L  =  1  jammed  BER  are  included  in  the 
next  section,  and  are  computed  as  a  special  case  of  the  L  >  1  BER  expression 
developed  there. 


138 


J.  S.  LEE  ASSOCIATES,  INC. 


4.2  ANALYSIS  FOR  DIVERSITY  SUM 

The  differential  detector  BER  performance  for  FH/CPFSK  using  diversity 
may  be  expressed  as 


PL(e)  =  \  •:  PL(e!  Ill)  +  2PL(e  Oil)  +  PL(el010)}  , 


where  L  is  the  number  of  hops  per  bit,  and 


L 

PL(e  xly)  =  Y.  (;)  (l-v)L"S-£  PL(e|xly,0 


In  this  expression,  the  possible  partial-band  jamming  events  are  indexed  by 
* ,  the  number  of  hops  jammed  out  of  L  for  a  particular  jamming  event.  As  a 
function  of  the  jamming  bandwidth  fraction  v,  the  probability  of  the  event 


Pr  ■  hops  jammed:  -  (I-*)*"' 


4.2.1  Derivation  of  Conditional  Error  Probability 

The  decision  statistic  for  L  hops/bit  diversity  is  the  sum  of  samples 
of  the  differential  detector  output.  As  was  shown  in  Section  4.1,  each 
sample  is  equivalent  to  the  difference  of  two  equally-scaled  noncentral 
chi-squared  random  variables,  when  the  cross-quadrature  correlation 
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j». 


coefficient  ■  is  zero: 


:k  —  ci[x2(2;dil)  -  x2{2;d.2)]  , 


(4.2-4a] 


where  the  jamming  condition  is  denoted  by 


0,  hop  unjammed 
1,  hop  jammed. 


(4.2-4b) 


From  Section  4.1,  we  have 


c0  4  "N  1_r2  ’  C1  =  Kc0; 


(4.2-5a) 


with  the  definition 


K  *  -- 


(Eb/No) 


(4.2- 5b ) 


Since  the  probability  of  error  is  not  affected  by  uniform  scaling  of  the 
samples,  we  may  simplify  matters  somewhat  by  using 


Cq  1,  c1  K. 


(4.2-5c) 


In  (4.2-4a),  the  noncentral i ty  parameters  d^.  and  d2 •  for  a  particular 
data  pattern  are  given  by 
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UM  -  r  Ww  cosA$  +  v7 1 - r 2  Ww  si n 

“01,02  ■  J - ----I - - -  ‘4'2-6a> 

and 

dll,12  =  ^d0rd02^/K’  (4.2-6b) 

with  U  and  W  being  SNR  parameters  and  being  the  (distorted)  differential 
phase  in  the  absence  of  noise.  The  SNR  is  related  to  the  bit  energy-to-noise 
density  ratio  by 


N 


LD 


(4.2-7) 


Since  chi-squared  variables  combine  to  form  chi-squared  variables  with 
higher  degrees  of  freedom,  the  diversity  sum  for  £.  hops  jammed  is  distributed 
as 


L 

2  =  zk~X:  [2(1-0  ;U-‘.)d01]  -  O[2(L-O:(L-Od02] 

k=l 

+  K -:(2.  ;idn)  -  KX2(2t  ^d12).  (4.2-8) 


4. 2. 1.1  Characteristic  Function  for  Sum 

The  characteristic  function  for  a  chi-squared  random  variable  with  2n 
degrees  of  freedom  and  noncentrality  parameter  d  is 


( v;n,d)  =  E  {  e' 


)  = 


_ 1 _ 

( l-2jv)n 


exp 


jvd 
1-2  jv 


(4.2-9) 
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Therefore,  the  characteristic  function  for  2  is 


*2(v)  =  *  .2[v;L-«.,(L-a)d0i]^  2[-v;L-fi.,(L-fi.)d02] 


x  v>  2[Kv;£,£d11]<f x2[-Kv;£,fcd12] 


=  [1+4v2]'(L_£)  [1+4K2v2]' 


■  exp 


.-DA  +  2,jv(  L-£] 
1  +  4v2 


-4Kv2£A  +  2jvRB 
1  +  4K2v2 


ft 

(4.2-10) 


in  which 


A  =  d0i  +  do2  =  2(UN  '  rWN  C0SA*)/(l-r’2) 


B  =  d01  '  d02  =  2WN 


(4.2-lla) 


(4.2-llb) 


4.2. 1.2  BER  From  Characteristic  Function 

The  cumulative  probability  distribution  for  a  random  variable  may  be 
written  in  terms  of  its  characteristic  function  as  follows  [20]: 


Pr '2  <  Z)  =  i  -  !  —  Im{^,(v)e‘jvZ  }  . 

2  7  v  z 


(4.2-12)  :>■: 


Therefore  the  BER  is  given  by 


P(e)  =  Pr{2  <  0)  =  \  ±  ^  lm{^z(v)} 


(4.2-13) 


-  r  _#»  j  *  j*  *_j»  - 


Application  to  the  conditional  probability  of  error  leads  to 


PL(e  jxly.t) 


1 

2 


[1+4v2]“(L'^  [1+4K2v2]'* 


-4v2(L-£)A  +  -4v2K*A 
l+4v2  1+4K2V2 


2v(L-fc)B  +  2v^B 
l+4v2  l+4K2v2 


The  integral  in  (4.2-14)  may  be  converted  to  one  with  a  finite  interval  of 
integration  by  using  the  following  change  of  variable: 


1  COSr 

2  1  +  sin:-  ’ 


for  which 

j  .  1  de 

'  '  2  1  +  sine  ' 

The  resulting  expression  is 

_  j  ^ 

PL(e,xiy,0  *|.i  [  ^  [1+(K2-1) (l-sine)/2]”C 

J7-/2 


exp 


-  j  (1-sine)  £  L-t 


2tK 


2+( K2-l ) (1-sine) 


B 

2 


cose 


.2-14) 


.2-15a) 


. 2- 1 5b ) 


x  sin 


_ 2t 

2+(K2-l)(l-sine) 


(4.2-16) 
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4.2.2  Numerical  Results  for  Differential  Detection 

The  BER  for  the  differential  detection  of  L  hops/bit  FH/CPFSK  in  partial- 
band  noise  jamming  was  computed  as 

PL(e)  =  \  (  PL(e| 111)  +  2PL(e|011)  +  PL(e|010)},  (4.2-17) 

where  . 

Pj_(e I  B.)  =  ^0)d-Y)L"\£  PL(elB,0  (4.2-18) 

*=o'  ' 

and  the  conditional  probabilities  P^(e|B_,t)  are  given  by  (4.2-16),  using 

A  =  2(U^  -  rWNcosAj)/(l  -  r2)  ( 4 . 2- 19a ) 

B  =  2WNs ini :/»■!“  F-'  .  (4.2-19b) 

The  subscript  "N"  denotes  that  U  and  W  are  calculated  using 

CNR  =  Eb/N0-  (4.2-19C) 

The  value  of  the  parameter  K  in  (4.2-16)  is 

K  =  (Eb/NQ)/(Eb/NT).  (4.2-20) 

For  the  results  shown  below,  we  have  used  Eb/N0  =  15  dB  and  D  =  WjpT  =  1. 
Figures  4.2-1  through  4.2-3  show  the  L  =  1  performance  of  the  differential 
detector  as  a  function  of  Eb/Nj,  and  parametric  in  y  ,  the  fraction  of  the  hop 
band  which  is  jammed.  The  figures  differ  in  that  h  values  of  0.70,  0.65,  and 
0.60  are  used,  respectively.  This  assortment  of  values  for  the  modulation 
index  was  used  because,  as  we  have  shown,  h  =  0.7  (and  D  =  1)  are  considered 
best  values  for  no  jamming,  while  in  [11]  h  =  0.6  (and  D  =  0.75)  are  said  to 
be  best  values  for  the  jamming  case. 

In  comparison  with  Figure  3.1-5,  we  first  note  that  all  three  of  the 
differential  detection  results  for  L  =  1  exhibit  a  higher  jammed  BER  than 
that  obtained  using  discriminator  detection.  This  is  not  surprising,  since 
differential  detection  is  known  to  be  less  effective  than  discriminator 
detection  without  jamming.  The  difference  in  performance  for  worst-case 


BIT  ENERGY  TO  JAMMING  DENSITY  RATIO  CdB) 

FIGURE  4.2-1  FH/CPFSK  DIFFERENTIAL  DETECTION  PERFORMANCE  IN 
PARTIAL-BAND  NOISE  JAMMING  FOR  L  =  1  AND  h  =  C 
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BIT  ENERGY  TO  JAMMING  DENSITY  RATIO  (dB) 


FIGURE  4.2-3  FH/CPFSK  DIFFERENTIAL  DETECTION  PERFORMANCE  IN 
PARTIAL-BAND  NOISE  JAMMING  FOR  L  =  1  AND  h  =  0. 
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jamming  is  about  0.75  dB,  for  the  parameter  values  used. 


Observing  now  the  differences  among  Figures  4.2-1,  4.2-2,  and  4.2-3,  we 


note  that  for  h  =  0.7,  there  is  a  definite  "curling  up"  of  the  parametric 


error  probability  curves  at  the  BER  value  of  10  .  This  indicates  that  the 


-6  -5 

noise-only  error  rate  is  between  10  and  10  for  h  =  0.7.  As  h  is  decreased 


to  0.65  or  0.6,  the  "curl"  straightens  out,  indicating  better  performances  in 


noise  for  these  values.  Thus  for  noise  only,  the  BER  for  the  differential 


detector  is  quite  sensitive  to  the  value  of  h.  For  example,  for  y  =  1,  a 


10  BER  is  obtained  for  the  following  total  SNR  values: 


Eb/Nr 


14.7  dB 


13.3  dB 


12.6  dB. 


This  2  dB  spread  in  wi de-band  noise  performance  is  not  echoed  in  the  worst- 


case  jamming  performance,  however.  The  Pur  relations  shown  in  each  figure 


reflect  only  a  0.54  dB  spread  in  performance  for  the  h  values.  Therefore, 


the  value  of  h  is  not  as  critical  in  worst-case  jamming  as  it  is  in  Gaussian 


noise. 


The  jammed  BER  performances  for  diversity  cases  of  L  =  2  and  L  =  3 


differential  detector  samples  added  are  shown  for  h  =  0.7  in  Figures  4.2-4 


and  4.2-5,  respectively.  The  worst-case  L  =  2  performance  is  0.36/0.275  = 


1.2  dB  worse  than  that,  for  L  =1;  for  L  =  3,  the  worst-case  performance  is 


0.50/0.275  =  2.6  dB  worse.  More  significantly,  these  worst-case  BER  results 


are  uniformly  worse  for  variation  in  E^/Nj,  having  the  same  inverse-linear 
dependence  upon  E^/Nj  as  the  L  =  1  case.  Thus  we  have  demonstrated  that 


simple  diversity  summing  of  differentially-detected  FH/CPFSK  samples  does 


not  yield  a  diversity  gain. 


From  [11]  we  know  that  diversity  combining  of  hard  decisions  does 
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FIGURE  4  2-4  DIVERSITY  SUM  FH/CPFSK  PERFORMANCE  IN  PARTIAL-BAND 
NOISE  JAMMING  FOR  DIFFERENTIAL  DETECTION  AND  L  =  2. 
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h  -  0.70 
D  =  1.00 
Ek/Nh  =  15.00  dB 
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)0  50. 00 

'BIT  ENERGY  TO  JAMMING  DENSITY  RATIO  (dB) 

FIGURE  4.2-5  DIVERSITY  SUM  FH/CPFSK  PERFORMANCE  IN  PARTIAL-BAND 
NOISE  JAMMING  FOR  DIFFERENTIAL  DETECTION  AND  L  =  3. 
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produce  a  diversity  gain  for  high  E^/Nq.  Therefore,  it  is  likely  that 
forms  of  soft-decision  combining  exist  which  are  better  than  hard-decision 
combining.  We  now  know  that  summing  of  samples  (soft  decisions)  is  not 
one  of  them. 
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APPENDIX  A 

DERIVATION  OF  M0DUL0-2^  DIFFERENTIAL 
PHASE  PDF  BY  DIRECT  METHOD 


From  Section  2.1.5  of  the  text,  the  joint  pdf  of  phases  is 
P,  (i ;  j^)  -  j  dRj  J  dR£  P^p(  R  i »  »R2»  ^2) 


where 


RiR2  -1 

PrP(R:.i:.R:,«:)  =  - — -  exp  Q. 

EP  4^(1--)  l--£  1 


2  ? 

r;+rp 

Qr  + 


URiR2  C0S^  -C2  +  ,0 


-  XR^  cos(:;  -v)  -  YR^  COS(t:  -w) 


-  2'-*':  ;  cos( :  ') . 


The  parameters  X,  v,  Y,  and  w  are  defined  in  (2.1-37)  of  the  text. 
Note  that  the  convention  we  have  adopted  is  that  subscript  "1"  refers  to 
quantities  at  time  t^  =  t$  and  "2"  refers  to  those  at  t£  =  t$  -  T.  Thus 
for  example,  1;  -  1  -4..  This  convention  is  opposite  to  that  in  [8j. 


.  v  V 
V  V 

•* 

N 

.  .  r» 
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A. 2  Integration  Over  Unneeded  Variable 

We  define  the  differential  phase  using  the  transformation  of 
variables 


X  =  Cl  -  $ 2 

y  = 


$2  =  x  +  y 

=  y 


with  unit  Jacobian.  Integration  of  y  over  a  2r  interval  yields 


P.,(*)  =  % 


where 


b(x)  =  -L 


di  cos^  du  u^  e  Au  Ig[u  •  b(x)] 
Jr 


— —  -\/  X2  +  Y2  +  2XY  cose  cos( x-v+w)  +  (Y2-X2)  sina 


.  2  . 1-. 


A. 3  Solution  for  Inner  Integral 


Using  the  integral 


(A-7) 


(A-8) 


(A-9) 


du  u3  e'Au  IQ(bu)  =  2  /  dy  y  e'2Ay  I^bv^) 


— ? —  p  ( 2  ■  1  •  — ) 

?  lr1 a a ) 

(2A;  1  1 


exp  (4bA  )  f1  +  $a]  » 


( A- 10 ) 
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we  obtain  the  pdf  expression 


»«<»>  ' 


da  -£252  exp{-C  +  ||  J  [1  +  ^  ]. 

(2A)i  HR  4A 


(A-ll) 


A. 4  Identification  of  Parameters 
From  before, 


2A  =  1  -  cosa  cos(x+;) 

=  1  -  (r  cosx  -  >sinx)cosa 


(A- 12) 


C  -  — 
L  4A 


=  — —  I  ci  +  ci  -  2u»/:io2  cos(Ac +:) 
l--:  1 


I  X2  +  +  2XY  cosa  cos(x-v+w)  +  (Y2-X2)sina 

4  1  -  (r  cosx  -  >sinx)  cosa 


U  +  V  sin  >  -  W  cosa  cos(x- ' 
1  -  (rcosx  -  >sinx)coSu 


( A- 1 3 ) 


l+fe=1-E  +  C 


=  1  -  E  +  2 


U  -  W(r  cos-^c  -  asintft) 
1  -  p2 


(A-I4) 


Therefore 


cosa  e 


SAd>-  XSinA 


PA,. V*; 
AC 


4  v 


u 


[l-(rcosx-  Xsinx)cosa] 


2  l 


l-L-r£. 


1  - 


Finally,  changing  the  sign  of  the  integration  variable  a,  and  that  of  x,  gives 
the  same  expression  as  (2.1-55).  The  first  sign  change  is  arbitrary.  The 
second  reflects  a  different  defining  convention  for x  . 
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APPENDIX  B 

A  GENERAL  PURPOSE  PROGRAM  FOR  THE  CPFSK  ERROR  PROBABILITY 

The  FORTRAN -77  program  listed  below  computes  the  CPFSK  bit  error 
probability,  including  clicks,  under  the  following  assumptions: 

(a)  Given  post-I.F.  SNR  (CNR),  modulation  index  h  =  2fdT, 
and  filter  bandwidth-time  product  D  =  W ^ pT ;  Gaussian- 
shaped  spectrum. 

(b)  Only  intersymbol  interference  effects  due  to  adjacent 
bits  are  significant;  eye  pattern  (differential  phase) 
components  contain  only  harmonics  with  frequencies 

*  Rb  =  l/T- 

In  addition  to  additive  white  Gaussian  noise,  the  receiver  is  assumed  to 
be  subject  to  Gaussian  noise  interference  with  Gaussian-shaped  spectrum  with 
given  SIR,  bandwidth,  and  center  frequency  defined  prior  to  the  receiver  filter. 
The  bandwidth  of  this  interference  is  specified  by  the  input  parameter 

£  =  Wj/Wj,-,  (B-: 


and  its  frequency  location  relative  to  the  signal  carrier  is  specified  by  the 
parameter 
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The  program  computes  the  error  using  the  F  function  given  by  (2.1-56),  and 
exploits  symmetries  to  calculate 

P(e)  =  \  (PQ  +  2Pj  +  P2),  (B- 

where 

Pi  =  \  [P(ex.Oy)  +  P(e|Xily.)] 

=  Pci[F(0;Al.)  -  \  F(M.-77;Acf.) 

+  |  F (r-Ai.;  -  Ai.)]  +  1-Pci;  (B- 

and  pci  =  exp{- jNci  | }  .  (B- 

The  patterns  are  specified  by 

x0  =  ^0  =  1 

*1  =  °,  yl  =  1  (B- 

x2  -  y 2  ~  o » 

and  the  pattern-dependent  parameters  are  listed  in  Tables  2.1-1  and  2.1-3  of 
the  text. 

For  no  jamming  or  partial-band  jamming,  we  use  the  input  parameters 
SIR  >>  1,  f  =  0,  and  6  >>  1,  with  the  SNR  in  the  program  being  Eb/NgD  when 
not  jammed,  or  E^/NTD  when  jammed. 
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